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ABSTRACT 


Three  boundary  value  problems  involving  discontinuities  in  a  parallel 
Pi  ate  waveguide  are  solved.  The  three  kinds  of  waveguide  discontinuities 
studied  are:  1)  the  metallic  step  discontinui ty,  2)  the  inhomogeneous  E-plane 
bifurcation,  and  3)  the  trifurcation.  The  generalized  scattering  matrix  tech¬ 
nique,  introduced  by  Mittra  and  Pace,  is  applied  to  solve  these  problems  .  The 
solutions  obtained  are  formally  exact  and  in  series  form.  No  restriction  is 
made  on  the  operating  f requel  cy  of  the  waveguide.  The  solutions  are  equally 
valid  for  oversized  waveguides,  in  which  several  modes  are  allowed  to  propagate 
down  the  guide  at  once.  The  results  of  numerical  computations  of  the  reflection 
coefficient  for  the  dominant  mode  in  the  waveguide  are  reported.  T.'ie  compu¬ 
tations  were  performed  on  a  digital  computer.  Where  possible,  the  author's 
results  are  compared  with  those  published  by  Marcuvitz,  Williams,  and  Cronson . 
Close  agreement  between  the  author's  results  and  those  of  the  above  authors  is 


noted . 
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1 .  INTRODUCTION 


Relatively  few  boundary  value  problems  with  applications  in  microwave  or 
antenna  engineering  can  be  solved  exactly.  Broadly  speaking,  boundary  value 
problems  which  are  amenable  to  exact  solution  fall  into  one  of  two  groups. 

With  the  first  of  these,  the  boundary  conditions  conform  to  one  of  the  co¬ 
ordinate  systems  in  which  the  scalar  Helmholtz  equation  is  separable.  In 
this  case,  the  partial  differential  equation  is  reduced  to  a  set  of  ordinary 
differential  equations,  the  solution  of  which  is  usually  quite  straightfor¬ 
ward.  A  sizable  volume  of  literature  concerning  the  separation  of  variables 
technique  is  available.  Morse  and  Feshbach*  provides  a  comprehensive  treat¬ 
ment  of  the  topic. 

There  exists  a  second  group  of  problems  which  can  be  solved  exactly  by 
means  of  integral  transforms  such  as  the  Fourier  and  two-sided  Laplace  trans¬ 
forms.  Often,  problems  of  this  sort  may  be  formulated  as  an  integral  equation 
of  the  Wiener-Hopf  type,  or  alternately  in  terms  of  certain  special  systems 
of  infinite  order  linear  algebraic  equations.  A  Wiener-Hopf  Integral  equation 
can  be  solved  by  the  application  of  the  Fourier  transform  and  certain  function 
theoretic  techniques.  This  method  for  solving  a  Wiener-Hopf  integral  equation 

is  called  the  'Wiener-Hopf  technique'.  Many  papers  on  the  application  of  the 

2 

Wiener-Hopf  technique  are  available  in  the  literature.  Noble  has  published 
an  excellent  text  concerned  with  both  the  theory  and  application  of  the  tech¬ 
nique.  The  exact  solution  of  systems  of  infinite  order  linear  algebraic 

3  4 

equations  by  function-theoretic  methods  is  discussed  by  Brillouin  $  Whitehead 

5  6  7 

Agronovich  et0  al  .  ,  Adonina  et  al.  ,  and  Hurd  and  Gruenberg  . 

Generally,  however,  the  solution  of  a  problem  can  only  be  formulated  in 
terms  of  a  differential  or  integral  equation,  or  system  of  equations,  which 


2 


can  be  solved  only  by  approximate  methods  Approximate  methods  are  many  and 

varied  They  include  variational  and  perturbational  technique s,  as  well  as 

f inite-dif ieience  methods  and  the  various  iteration  procedures  commonly  used 

to  solve  integral  equations  .  Again;  the  volume  of  literature  on  the  subject 

8  9 

is  enormous .  Hartree  and  Householder  have  published  well-kn  vn  texts  on 

10 

numerical  analysis,  Goertzel  and  Tralli  is  representative  of  the  general 
references  available  on  mathematical  physics,  A  wide  range  of  topics  is 
covered>  including  chapters  on  perturbation  of  eigenvalues^  variational  esti- 
mates>  etc, 

A  new  technique  for  solving  a  class  of  boundary  value  problems  is  dis¬ 
cussed  in  this  thesis.  It  will  be  referred  to  hereafter  as  the  generalized 
scattering  matrix  technique  for  reasons  which  will  be  made  clear  in  the  sub- 
sequen*  discussion.  The  application  of  this  technique  makes  it  possible  to 
derive  a  formally  exact  solution,  in  series  form,  to  problems  for  which  only 
approximate  solutions  have  been  possible  before . 

It  is  believed  that  the  generalized  scattering  matrix  technique  should 
have  a  broad  range  of  applicability-  The  purpose  of  this  thesis,  however^ 
is  to  demonstrate  its  usefulness  for  solving  certain  boundary  value  problems 
associated  with  discontinuities  in  a  parallel  plate  waveguide  Specifically, 
three  distinct  problems  are  discussed.  They  are  the  E-plane  metallic  step 
discontinuity  (Figure  1)7  the  inhomogeneous  E-plane  bifurcation  (Figure  2) } 
and  the  trifurcated  waveguide  (Figure  3>,  The  inhomogeneous  E-plane  bifur¬ 
cation  is  an  ordinary  bifurcated  waveguide  modified  by  placing  a  dielectric 
in  one  of  the  smaller  ducts  of  the  waveguide,  i,e«,  with  reference  to  Figure  2, 

the  dielectric  is  placed  in  region  B 

.  1  1  r  ,  12  13  .14  .  u  15 

I.ewm  ,  Collin  ,  Ghose  ,  Durrani  ,  and  Harvey  provide  a  survey  of 


waveguide  theory  and  existing  techniques  for  solving  waveguide  discontinuity 
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Figure  1.  The  S-plane  metallic  step  discontinuity. 


REGION  A 


REGION  B 


REGION  C 


Figure  2.  The  inhomogeneous  E-plane  bifurcat 


REGION  B' 


REGION  A  I 

j  REGION  B" 

1 - 

!  REGION  C 


Figure  3.  The  trifurcated  waveguide 
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problems.  Some  papers  of  note  specifically  concerned  with  the  three  wave¬ 
guide  discontinuities  discussed  in  this  thesis  have  been  published.  They  are 
discussed  in  Chapters  4,  5 f  and  6  in  which  the  problems  themselves  are  dis¬ 
cussed  . 

This  concludes  the  introduction.  In  the  next  chapter.,  the  generalized 
scattering  matrix  technique  is  discussed. 
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2.  THE  GENERALIZED  SCATTERING  MATRIX  TECHNIQUE 

16 

In  this  section  of  the  paper,  a  new  technique  for  the  solution  of  a 
class  of  boundary  value  problems  arising  in  electromagnetic  theory  is  pre¬ 
sented.  Although  the  technique  should  be  applicable  to  other  kinds  of  prob¬ 
lems,  too,  the  generalized  scattering  matrix  technique  is  explained  here  by 
relating  it  to  problems  involving  a  class  of  waveguide  discontinuities.  In 
particular,  the  technique  will  be  applied  to  three  boundary  value  problems 
involving  a  parallel  plate  waveguide  configuration.  They  are  the  inhomogene¬ 
ous  E-plane  bifurcation,  the  metallic  step  discontinuity,  and  the  trifurcated 
waveguide.  Only  brief  mentio.i  of  these  problems  is  made  in  this  section, 
however.  The  detailed  solutions  follow  in  later  sections. 

The  key  to  the  technique  is  the  identification  of  an  auxiliary  problem 
associated  with  the  particular  problem  to  be  solved.  Each  of  the  three  prob¬ 
lems  discussed  in  detail  in  this  thesis  has  a  common  auxiliary  problem.  This 
is  the  boundary  value  problem  associated  with  a  semi-infinite  bifurcation  in 
a  parallel  plate  waveguide  (refer  to  Figure  4).  The  geometry  of  the  auxil¬ 
iary  problem  is  such  that  it  can  be  modified  in  a  straightforward  manner  so 
as  to  be  made  identical  with  the  geometry  of  the  original  problem.  For  in¬ 
stance,  if  region  B  of  the  bifurcation  is  filled  with  dielectric,  then  the 
resulting  configuration  is  just  that  of  the  inhomogeneous  E-plane  bifurcated 
waveguide.  The  dielectric  can  be  regarded  as  a  termination  or  load  placed  in 
region  B.  Likewise,  the  appropriate  modification  in  the  case  of  the  metallic 
step  discontinuity  is  made  by  placing  a  perfectly  reflecting  wall  in  region  B 
flush  with  the  plane  of  the  junction  (z  -  0).  In  the  case  of  the  trifurca¬ 
tion,  the  auxiliary  problem  is  modified  by  placing  a  second  bifurcation  in 


region  B, 


figure  4.  The  inhomogeneous  E-plane  bifurcation 


The  auxiliary  problem. 


(Use  this  Figure  4  in  place  of  the  page  8  Figure  4  in  text) 


BLANK  PAGE 
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The  auxiliary  problem  is  characterized  by  the  fact  that  it  can  be  solved 

exactly.  The  problem  of  the  semi-infinite  bifurcation  can  be  solved  exactly 

by  any  one  of  several  methods.  For  instance,  it  can  be  handled  by  the 
17  7 

Wiener-Hopf  or  function-theoretic  technique  ,  or  the  direct  solution  of  an 

1 8 

infinite  set  of  linear  algebraic  equations 

The  way  in  which  the  auxiliary  problem  is  used  to  solve  the  three  above 
mentioned  problems  will  be  discussed  next. 

2ol  Development  of  the  Series  Expansion 

It  is  common  practice  to  think  of  the  bifurcated  waveguide  in  terms  of  a 
transmission  line  analogy.  If  only  the  dominant  mode  of  the  guide  is  allowed 
to  propagate,  as  usually  is  the  case,  the  bifurcated  waveguide  is  regarded  as 
a  3-port  network.  One  port  of  the  network  is  associated  with  the  propagating 
mode  in  each  of  the  three  regions  A,  B,  and  C,  A  port  may  be  regarded  as  the 
terminals  of  a  transmission  line.  If  the  problem  of  multi-mode  propagation  is 
under  consideration,  then  the  network  is  regarded  as  an  V-p0rt  system,  N 
being  the  total  of  thfv  propagating  modes  in  +he  three  regions 

In  the  neighborhood  of  the  edge  of  the  discontinui! y,  an  infinite  number 
of  evanescent  modes  are  excited  by  the  diffraction  of  an  incident  plane  wave. 
The.  effect  of  these  evanescent  modes  can  be  represented  in  terms  of  lumped 
reactances  in  the  equivalent  circuit  4s  these  modes  in  physical  terms  repre¬ 
sent  stored  energy.  Of  course,  the  numerical  values  of  these  lumped  re¬ 
actances  must  be  determined  by  solving  the  boundary  value  problem. 

The  network  of  lumped  reactances  and  transmission  lines  can  be  concisely- 
described  in  mathematical  terms  by  means  cf  a”  impedance,  admittance,  or 
scattering  matrix.  The  order  of  any  of  these  matrices  will  be  N  where  N  is 
the  total  number  of  propagating  modes  in  the  three  regions  A,  B,  and  C 
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In  the  technique  introduced  in  the  thesis  the  bifurcated  guide  is  re¬ 
garded  as  a  3N-po»'t  •'etwork^  v  bei*-g  i-rinitc-lv  large,  A  port  of  *he  network 
is  assigned  to  each  of  the  propagating  and  evanescent  modes  If  the  ports  of 
the  network  are  termi-ated  with  the  proper  set  of  ref  le  ct  *nces,  a  network 
corresponding  to  either  the  step  discontinui Xy j  the  inhomogeneous  bifurcated 
waveguide,  or  the  trifurcation  is  achieved, 

Hie  concept  of  a  scattering  matrix  of  infinite  order  is  introduced « 

While  the  mechanism  of  its  application  is  convent ional,  it  differs  from  the 
scattering  matrices  ordinarily  cefined  in  the  literature,  For  one  thing;  the 
concept  of  the  scattering  coefficient  is  extended  to  cover  evanescent  modes. 
The  following  discussion  is  concerned  with  the  derivation  of  the  generalized 
scattering  matrix  of  infinite  order  as  applied  to  waveguide  discontinuity 

problems  of  the  kind  discussed  ir  this  thesis 

cm 

With  reference  to  Figure  4,  let  S  a  ■  A,  P.  or  C  represent  the  self¬ 
scattering  matrices  of  the  auxiliary  problem,  Let  Sa^  a  A  B  or  C  and 
p  A,  B,  or  C  but  a  /  (3^  be  the  mutual-scattering  matrices. 

The  interpretation  of  the  scattering  matrices  is  as  follows.  Consider 
that  regions  B  and  C  are  terminated  i"  reflectionless  loads.  Then  if  the 
n**1  transverse  magnetic  mode  is  incident  ir*  region  A  fields  will  be  reflected 
in  region  A  and  transmitted  to  regions  B  and  C,  These  fields  can  be  deter¬ 
mined  by  solving  the  auxiliary  problem.  The  resultant  electromagnetic  fields 
are  expressible  entirely  in  terms  of  transverse  magnetic  modes.  Thus,  the 
total  electric  field  is  expressible  in  terms  of  the  total  H  field,  Ir  turn, 

y 

H  can  be  written  in  *erms  of  eigenfurct  io^  expansions  with  constant  co- 

y 

efficients  appropriate  to  regions  A,  P  a^d  C,  If  *his  thesis/  the  com¬ 
ponent  of  the  transverse  magnetic  field  is  expanded  in  each  of  the  three 
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regions  A  B  and  f  ir  terms  of  Fourier  cosine  series  of  the  form  £  d  * 

n  n^n 

,  .  nx  _  ^  r ( x -  a ' 

where  in  rogion  A  d>  cos  -  in  region  B  d>  cos - - -  and  m 

B  '  vn  a  b 


_  ,  ^  nx 

region  C ,  d>  -■  cos  - 


The  mode  coefficient  of  rhe  m  mode  referred  to 


the  plane  of  the  junction  (z  0)y  appearing  in  the  expansion  of  will 

be  defined  as  the  ’amplitude'  of  the  mode.  If  the  amplitude  of  the  n  mode 

incident  in  region  A  is  one,  or  in  other  words,  the  incident  field  is 
If  nx 

given  by  cos  - ,  the  total  H  field  reflected  in  region  A  at  z  0  is  given 

a  *  y 

AA  ft  mx 

by  -  £  S  cos  -  .  The  total  H  field  transmitted  say  to  region  B  at 

m -  0  mn  a  y 

^  BA  Ttm(x-a)  „  th 

z  0  is  given  by  zL  S  cos  - - - .  hus,  if  the  amplitude  of  the  n 

m  0  mn  b  J 


mode  incident  in  region  A  is  one(  the  amplitude  of  the  m  mode  scattered  in 

AA  b 

region  A  will  be  -S  and  the  amplitude  of  the  m  mode  transmitted  to  region 

mn 

BA  AA  BA  AA  BA 

B  will  be  S  »  S  and  S  are  the  general  matrix  elements  of  S  and  S  , 
mn  mn  mn 

respectively «  The  other  matrices  are  defined  in  a  similar  manner 

Note  the  scattering  coefficients  are  defined  in  this  thesis  in  a  manner 
consistent  with  the  sign  convention  followed  in  electromagnetic  theory  when 
defining  the  reflection  and  transmission  coefficients  in  terms  of  the  voltage 
or  transverse  electric  field  rather  than  i"  terms  of  the  current  or  trans¬ 
verse  magnetic  field.  Since  the  coefficients  of  the  eigenfunction  expansion 
of  H  are  used  directly  to  define  the  various  scattering  coefficients  care 

y  J  ■’ 

must  be  taken  to  assign  the  proper  sign  to  the  ratio  of  the  amplitude  of  the 
scattered  mode  to  the  amplitude  of  the  incident  mode  if  the  definitions  are 


to  be  consistent  with  this  sign  convention  Thus,  the  amplitude  of  the  m 

AA  AA 

mode  comprising  H  m  region  A  is  given  by  -S  and  not  by  just  S 

y  mn  mn 

It  should  be  noted  that  ordinarily  the  mode  amplitudes  are  normalized  so 


that  a  propagating  mode  carries  unit  power.  Hcwever/  since  the  scattering 
matrix  has  been  generalized  to  include  evanescent  modes,  it  is  inappropriate 


12 


to  normalize  the  mode  amplitudes  in  the  usual  magnet  7"  ♦his  thesis  the 
scattering  coefficients  are  defined  a-  the  ratio  cf  tie  amplitude  of  a  scat¬ 
tered  mode  tc  the  amplitude  of  tie  incident  mode  which  is  taken  to  be  unity. 
One  consequence  of  this  definition  however  is  t~*t  the  various  scattering 
matrices  are  non-symmet ric 

Suppose  now  that  region  B  is  modified  by  placi-g  an  obstacle  in  it  such 
as  a  perfectly  reflecting  wall  I*"  terms  of  the  ^~port  network  representa- 
tiony,  the  load  in  region  B  car  be  represented  by  a  reflectance  matrix  T  and 

D  0  0 

a  transmission  matrix  w  „  The  meaning  of  r  *na  4  can  be  interpreted  as 
follows.  Let  t  be  a  vector  associated  with  the  field  distribution  of  a  wave 
progressing  toward  positive  z  if  region  P.  The  e':emerts  of  the  vector  X~  of 
coarse ,  are  the  amplitudes  of  the  modes  of  the  eigenfunction  expansion  of  the 

arbitrary  field  in  region  P  The  ampli’udes  are  referred  to  +he  plane  z  0, 

B-  p.~- 

If  this  is  so,  then  r  t  gives  the  reflected  field  a"d  4  t  gives  the  trans- 

B  — 

mitted  fields  both  referred  to  the  pl*ne  of  the  junction  at  z  r  •  is 

progressing  in  region  B  *oward  f'egitive  z  away  from  the  load.  The.  elements 

0 _ 

of  r  t  are  the  amplitudes  of  the  trades  in  the  reflected  field.  Similarly, 
4^  is  progressing  in  region  P  toward  positive  z  i"d  *he  elements  of  4**T  are 
the  amplitudes  of  the  modes  comprising  tie  transmitted  field.  As  an  example^ 

if  a  perfectly  reflecting  wall  in  region  B  is  flu**-  with  *he  plane  of  the 

B  -■ 

junction,  T  -I  where  ’  is  the  identity  matrix  a"o  4  0 

0  B 

With  r  md  4  defineo  consider  the  situation  shown  n  Figure  5>  The 
usual  boundary  value  problem  associated  with  this  kind  of  geometry  involves 
a  monochromatic  signal  excited  in  some  region  of  *ne  guide  and  propagating 
toward  the  junction  O^c  is  interested  in  determining  » W  fields  diffracted 
by  the  junction  or  at  !e  i-t  some  part  of  them  F or  example,  the  reflection 
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REGION  A 


REGION  B 

load  rB 


REGION  C 


Figure  5.  Auxiliary  problem  modified  by  a  load 
placed  in  region  B. 
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and  transmission  coefficients  for  the  dominant  mode  are  often  sought.  One  way 
to  express  the  desired  fields  is  in  terms  of  their  eigenfunction  expansions 
with  constant  coefficients.  The  generalized  scattering  matrix  technique  en¬ 
ables  one  to  write  the  coefficients  of  the  eigenfunction  expansion  in  terms  of 
the  scattering  coefficients  of  the  auxiliary  problem  and  the  load  in  region  B. 
The  derivation  of  the  relationship  between  the  eigenfunction  coefficients  and 
the  scattering  matrices  follows . 

Suppose  region  A  is  excited  by  a  TEM  mode.  Suppose  also  that  region  B  is 

B 

terminated  by  a  load  characterized  by  the  reflectance  matrix  r  „  The  TEM  mode 

will  be  scattered  and  reflected  first  by  the  bifurcation  in  the  waveguide.  The 

field  reflected  into  region  A  can  be  characterized  by  the  vector  r  where 
—  AA  —  —  T 

rQ  =  S  a.  By  definition  a  --  (1,  0,  0;  .,.)  where  T  means  the  transpose. 

(Actually,  in  this  particular  case  rQ  0).  Furthermore  a  wave  will  be  trans- 

—  —  BA  — 

mitted  to  region  B.  Let.  t  characterize  this  wave  where  t  =  S  a.  The  wave 

o  o 

will  be  reflected  by  the  load  in  region  B.  The  reflected  wave  is  character- 
_B— 

ized  by  1  tQ.  The  reflected  wave  progresses  toward  the  junction  in  the  nega¬ 
tive  z  direction  where  it  is  diffracted  by  the  edge  of  the  bifurcation.  A 

field  r,  ,sABr*r  is  transmitted  to  region  A  and  a  field  t,  SBBr®  t  is 
1  o  1  o 

reflected  in  region  B.  This  field  will  also  be  reflected  back  by  the  termi¬ 
nation  and  this  process  of  multiple  reflection  will  be  continued.  All  of  the 
contributions  in  region  A  due  to  this  process  can  be  written  symbolically  as 


S 


in 


00 

-  20  7 

n-0  n 


AB_B  — 

s  r  t 


s^r8 


sBBrBT 


(i) 


This  is  recognized  as  a  Neumann  type  series.  The  Neumann  series  can  be  sum¬ 
med  in  the  usual  manner,  and  Equation  (1>  can  be  written  as 


S  -  S**  a 
in 


AB  ,-B  BB  ~  B  -1  BA— 

s  r°  (i  -  s  r  )  1  S  a 


<2) 
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—  I 

S.  is  given  by  (R  ,  -A,.  -A„,  . . .)  where  R.  is  the  voltage  reflection  co¬ 
in  A  12  A 

efficient  for  the  dominant  mode  and  the  A  *s  are  the  coefficients  of  the 

n 

Fourier  series  expansion  of  in  region  A .  Note  that  the  difference  in  the 

signs  between  R.  and  the  A  *s  above  is  due  to  the  aforementioned  definitions 
A  n 

of  the  scattering  coefficients.  The  convergence  of  the  Neumann  series  is 
discussed  in  Section  2.2. 

The  discussion  of  the  derivation  of  Equation  (2)  in  the  preceding  para¬ 
graph  depends  on  an  intuitive  understanding  of  the  physical  processes  involved 
in  the  chain  of  multiple  reflections.  In  order  to  provide  a  lucid  explanation^ 
it  was  tacitly  assumed  that  the  load  in  region  B  was  recessed  some  arbitrary 
distance  6  into  region  B  and  away  from  the  edge  of  the  bifurcation.  Refer  to 
Figure  6°  If  indeed  this  is  the  case;  then  the  proof  of  the  convergence  of 
Equation  (1)  is  very  much  simplified  since  all  of  the  higher  order  evanescent 
modes  would  be  damped  out  quite  rapidly.  However^  this  crutch  is  not  neces¬ 
sary.  In  the  limit 3  6  can  be  zero  and  it  will  be  shown  that  the  series  even 
then  is  convergent .  Indeed^  convergence  is  proven  without  reference  to  the 
specific  expressions  for  the  various  scattering  coefficients. 

Using  reasoning  similar  to  that  employed  in  deriving  Equations  (1)  and 
(2),  an  expression  for  the  fields  in  say  region  C  can  be  written  as 


-  CA  -  _CB  B  .  BB  „B.-1  BA  - 

sc.a  • s  a  +  s  r  "  - s  r  ’  s  a 


(3) 


where  S  -  (T  ,  C  ,  C  ,  ...)  .  T  is  the  transmission  coefficient  for  the 

dominant  mode  from  region  A  to  region  C  and  the  G^'s  are  the  higher  order  co¬ 
efficients  of  the  Fourier  series  expansion  for  the  solution  field  in  region  C 
In  a  similar  manner^  the  mode  coefficients  of  the  eigenfunction  expan¬ 


sion  of  Hy  in  region  B  can  be  written  as 

-  JB  /T  BB  ,-l  0BA  - 
SgA  -  <x>  (I  -  S  r  )  S  a 


(4) 
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REGION  A 


REGION  C 


z*0 
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—  T 

where  --  (T  ^  B ^  . .»)  .  is  the  transmission  coefficient  for  the 

dominant  mode  from  region  A  to  region  B  and  the  B_  ' 3  are.  the  higher  order 
mode  coefficients  of  the  Fourier  series  expansion  for  in  region  B. 

Equations  (2)#  (3) /  and  (4)  apply  to  the  specific  case  of  a  signal  ex¬ 
citing  region  A  and  a  load  placed  in  region  Bo  However,  similar  equations 
can  easily  be  derived  for  other  situations  *  Also^  note  that  loads  can  be 
placed  in  two  regions  in  the  same  problem .  For  instance^  the  problem  of  the 
capacitive  diaphragm  in  a  waveguide  can  be  solved  through  the  device  of  plac¬ 
ing  a  magnetic  wall  in  region  C  and  an  electric  wall  in  region  Bo  The  method 

19 

of  multiple  scattering  has  been  applied  to  this  problem  for  the  case  of  a 
semi-septum.  The  solution  was  shown  to  be  quite  st raightf orward » 

2«2  The  Proof  of  the  Convergence  of  the  Neumann  Series  Expansion 

In  this  section,  the  convergence  of  the  Neumann  series  expansion  given 


by 


,  BB  _B  -1  BB  BB  BB  ^ B 

(i  -  s  r  >  ...  i  ~  s  rD  +  s  r  s  r 


(5) 


is  proven.  First,  the  convergence  of  the  series 


,  BB  -1  BB  BB  BB 

(I  -  S  )  I  f  S  ;  S  S 


(6) 


is  demonstrated.  Next;,  it  is  demonstrated  that  if  series  16)  is  convergent > 
series  ( 5 i  is  also  convergent. 

Let  a  and  b  be  the  amplitudes  of  the  ntb  mode  incident  and  reflected 
n  n 

in  region  B^  defined  at  the  plane  of  the  junction  (z  0)  >  Furthermore,,  let 

t  h 

and  i^  be  the  voltage  and  current  for  the  n  mode  in  region  B,  defined  at 


1  _  (p)  |2  (pi 

—  Z  a  I  .  where  Z 

2  n  n|  >  n 


2u* 


For  | 

propagating  modes. 

w 

to 

1 

'Eh 
i b 

2 

o  Similarly, 

also  with  b  ~  0,  one  has  ~  e  i*--^Z 
n  *  2  n  n  2  n 


(e) 


2  .  _  (e) 

a  .  where  Z 
n  ‘  n 


2wt 


Al 


fffn  "  2 

T"  *k 
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For  Propagating  nodes,  eQ  and  i  can,  in  general,  be  expressed  in  terns  of 


a  and  b  as 

n  n 


e  *  Z  <P)  (a  +  b  ) 
n  n  n  n 


(7) 


and 


i  =  (»  “  b  ) 

n  n  n 


For  evanescent  nodes,  e  and  i  can  be  expressed  as 

n  n 


e  =  -JZ  <e)  (a  +  b  ) 
n  n  n  n 


(8) 


(9) 


and 


i  «  (a  -  b  ) 
n  n  n 


(10) 


With  the  reference  plane  chosen  in  region  B  at  the  plane  z  o,  regions 
A  and  C  can  be  grouped  together  and  regarded  as  a  termination.  A  relation¬ 
ship  between  the  terminal  voltages  and  currents,  the  stored  energy,  and  the 

20 

power  delivered  to  the  termination  is  derived  by  Jlontgomery,  et  al . 


eo 

1  £  e  i  * 

2  n=to  n  n 


2jw  (WH  -  Wg)  +  P 


(11) 


where  Wu •»  average  magnetic  energy  stored  in  the  te rain at ion 

W  =2  average  electric  energy  stored  in  the  termination 

•  \  £ 

p  is:  average  power  delivered  to  the  termination. 

To  simplify  the  following  discussion,  assume  that  only  the  TEM  mode 
propagates  in  region  B.  Substituting  (7),  (8),  (9)  and  (10)  into  Equation 
(11),  one  obtains 

•  <P> 


(a  t  b  )  (a  -  b  )  -  j  2L  Z 
o  o  o  o  n  >  l  i 


(e)  *  * 

(a  +  b  )  (a  -  b  ) 
i  n  n  n  n 


4jw  <WH  -  Wg)  *  2P 


(12) 
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Equating  the  real  parts  of  Equation  (12),  one  obtains 

2  (P>  (la  I2  -  lb  I2)  -  j  2  Z _<e>(a  *  b  -  a  b  *)  =  2P  (13) 
o  I  ol  |  o  |  n=l  n  n  n  n  n 

Equating  the  Imaginary  parts  of  Equation  (12).  one  derives  the  equation 


Zo<P>  (*obo*  -  *o\>  *  J  n?l  Z„(e)  <k|*  '  K  F>  “  «  <*E  -  WH>  (14> 


—  BB  —  —  BB 

Now  let  b  -  S  a  where  a  is  an  eigenvector  of  S  ,  so  that  if  \  is  an 


eigenvalue 


b  =  \  a 


(15) 


b  ~  \  a 
n  n 


(16) 


Substituting  Equation  (16)  into  Equations  (13)  and  (14).,  one  obtains 


Z  <P)  a  2  (1  -  \|2)  +  j  2  <\*  -  \)  Z  <e)  a  |2  ,  2P 

o  o  w  n^l  n  n 


(17) 


Z  (P)  a  2  (\*  -  \)  f  j  2  (1  -  \  2)  Z  (e)  a 

o  o  Jn~l  n  ] 


n|2  *  4  J“<*E  *  »H>  <18> 


Equations  (17)  and  (18)  can  be  treated  as  two  equations  in  the  two  un¬ 
knowns  (\  -  \)  and  (1  -  |\.|2)o  Solving  for  (1  -  |\|2),  one  gets 


2P  Z  (p)  la  I2  -  4W(W  -  W  )  2.  Z  (e)  |a  |2 

1  _  I  2 _ _ o  I  oj _ E  H  n  1  n  |  n| 

(Z  (p)  la  j2)  e  (  2  Z  (p)  la  j2) 
o  o  n  1  n  n 


(19) 


Note  that  Z  ^P^  and  Z  are  positive  real  quantities  and  >  W..  for 

o  n  EH 

the  case  of  Ttl  modes ,  Thus, 


\  -  1 


(20) 
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A  necessary  and  sufficient  condition  for  the  convergence  of  the  Neumann  series 

I  +  II  +  M8  +  ...  is  that  the  eigenvalues  of  li  satisfy  the  inequality  |\  |<  1. 

BB 

Th  the  convergence  of  (I  +  S  )  is  proved. 

BB  BB  BB 

The  proof  of  the  convergence  of  I  +  S  +  S  S  ...  for  the  case  of  TE 

•odes  is  formally  the  same  as  the  proof  given  above  for  TM  modes.  In  this 

th 

case  V  >  W  ,  but  the  voltage  e  for  the  n  non-propagating  voltage  is  nor- 
H  a  n 

(e)  (e) 

■a 11 zed  differently.  In  this  case,  e  =  J  Z  (a_  +  b  ),  Z  >0. 

#  n  n  n  n  n 

There  is  an  alternate  condi tion,  necessary  and  sufficient,  for  the  con- 

21  — 

vergence  of  a  Neumann  series  (see,  for  instance,  Friedman  ).  If  a'  is  an 

— ?  —  2 

arbitrary  vector  «nd  b*  =  Ha',  then  the  Neumann  series  I  +  If  +■  M  +  ...  con¬ 
verges  absolutely  if  jb'J  <  |a*|  where  |  a*|  is  finite.  Thus,  from  the  proof 

BB  BB  BB 

of  the  convergence  ofl+S  +  S  S  +  ...,  one  already  knows  that  if 
b1  =  SBB  a',  then  Jb*|  <  Ja*|  .  Now,  let  a*  =  T®?'.  Then,  if  b"  •=  SBBr^a", 


B  BB 

Now,  one  can  follow  the  same  argument  for  r  as  was  used  for  S  to  de¬ 
rive  Equation  (19),  again.  In  this  case,  however,  it  is  possible  that  P  -  0 

__B 

and  W_  =  W  as  would  be  the  case  if  I  *  *  I.  Hence,  if  *1  is  any  eigenvalue 
EH  ”  . 

of  rB,  |r| |  <  1,  and  it  follows  that  |r®  a’j  <  (a'l  .  Thus,  |b"|  <  j  a  J  , 

BB— B  BB-B  BB— B 

showing  that  the  Neumann  series  I  +  S  IT  +  S  r  S  *  +  . . .  is  absolutely 

convergent  and  can  be  summed  to  (I  -  sBBrV  k.  The  proof  is  essentially  the 
same  if  the  case  of  several  modes  propagating  is  considered. 

In  Chapters  4,  5,  and  6,  (I  -  S88!*8)  is  truncated  and  then  inverted.  The 
Neumann  series  expansion  was  used  in  computing  the  inverse.  For  all  cases 
considered,  it  was  found  that  the  Neumann  series  can  be  truncated  after  20  to 
30  terms.  The  computation  of  the  inverse  of  (I  -  S88!3)  using  the  Neumann 
series  expansion  is  straightforward  even  when  the  rank  of  the  truncated  matrix 
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BB  R 

is  relatively  large.  Of  course,  if  the  rank  of  the  truncated  matrix  (I  -  S  F  ) 

is  sufficiently  small,  then  the  matrix  can  be  inverted  directly 

2o3  General  Comments  on  the  Technique 

The  introduction  of  an  auxiliary  problem  of  the  kind  discussed  above  as 

an  aid  to  the  solution  of  certain  kinds  of  boundary  value  problems  is  rather 

new,  A  search  of  the  literature  has  uncovered  only  one  paper  in  which  the 

22 

author  uses  a  similar  device,  W>  E,  Williams  uses  *he  Laplace  transform  in 
the  formulation  of  the  step  discontinuity  problem.  He  applies  the  Wiener-Hopf 
technique  and  derives  an  associated  set  of  infinite  order  linear  algebraic 
equations.  The  auxiliary  problem  is  introduced  by  Williams  as  a  preliminary 
step  to  solving  the  infinite  set  of  equations.  This  is  to  be  contrasted  with 


the  use  of  the  auxiliary  problem  as  discussed  in  this  thesis,  Using  the 


generalized  scattering  matrix  technique,  the  solution  to  the  problem  is  ex¬ 
pressed  in  terms  of  the  scattering  coefficients  of  the  auxiliary  problem  as 
a  rapidly  convergent  series. 

Also,  the  applicability  of  the  generalized  scattering  matrix  technique 
should  be  rather  broad.  The  method  is  not  necessarily  restricted  to  wave¬ 
guide  problems.  For  example,  it  is  suggested  tha*  the  problem  of  the  dif¬ 
fraction  of  a  plane  wave  by  a  thick  conducting  half-plane  may  be  solved  by 
means  of  the  new  technique  The  auxiliary  problem  suggested  is  the  boundary 
value  problem  associated  with  a  pair  of  parallel,  semi -inf i m te  plates  in 
free  space  (Figure  7), 

A  brief  list  of  problems  suggested  for  future  s*udy  is  given  in  Chapter 


7  of  the  thesis „ 
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lb) 

Figure  7b.  The  proposed  auxiliary  problem:  parallel-plate 
waveguide  In  space. 
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3.  DERIVATION  OF  THE  SCATTERING  COEFFICIENTS 

The  elements  of  the  self-scattering  matrices  S40,  a  A  B,  or  C.  and  the 

a3 

mutual-scattering  matrices  S  ^  a.  A  c*  or  C  and  (3  A,  P  or  C  but  a  /  (3, 
are  determined  by  solving  the  boundary  value  problem  associated  with  a  semi¬ 
infinite  bifurcation  in  a  parallel  plate  waveguide.  The  problem  must  be 
solved  for  an  arbitrary  TM  mode  incident  from  one  of  the  three  regions  A, 

BJ  or  C. 

It  may  be  shown  for  the  problem  under  consideration  that  the  only  modes 

excited  by  the  discontinuity  with  an  arbitrary  TM  mode  incident  are  the  TM 

7  7  no  no 

modes o  The  non-vanishing  field  components  can  be  derived  from  a  single  scalar 
function  jS  (x;  z)  which  is  identical  to  the  y-component  of  the  H-field,  Using 
the  coordinate  system  illustrated  in  Figure  4  the  three  components  of  the 
field  can  be  written  as 

H  --  f)  (21 ' 

y 


:  -1- 
x  Oz 


(22) 


and 


-1 


Z  3x 


(23) 


-ju;t 

A  harmonic  time  variation  of  the  type  e  is  assumed  throughout  The  scalar 
f unct ion  *(x,z>  must  satisfy  the  two-dimensional  Helmholtz  equation 


d2  d2  2 

2  2  k 

ax  az 


4)  o,  k 


2TT 


(24) 


together  with 


e-«. 


x  0 .  a  a  1 1  z  and  x  -  t  z  0 


(25  > 
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and  the  edge  condition 


|v  <$>  |  0(d~1/2)'  d 

-p  ♦  (X  -  c)*p2 


(26) 


where  d  =|z  +  (*  "  C)J 

Initially,  it  is  assumed  that  the  waveguide  is  excited  from  each  of  the  three 
regions  A,  B,  and  C  by  an  arbitrary  transverse  magnetic  mode  of  order  I,  q, 
and  r,  respectively.  It  is  assumed  that  the  waveguide  dimensions  are  such 
that  each  of  the  incident  modes  is  a  propagating  mode.  Thus,  in  general, 
there  will  be  several  propagating  modes  in  each  of  the  three  regions  of  the 
waveguide.  Each  of  the  propagating  modes  must  satisfy  the  radiation  condition 
at  infinity. 

3.1  Derivation  of  the  Inf inite  Sets  of  Equations 


In  each  of  the  three  regions  labeled  A,  B,  and  C  the  function  ,  (x,z)  can 
be  written  in  terms  of  the  appropriate  eigenfunction  expansion,  in  this  case 
a  cosine  series .  Thus,  in  region  A 


-  A  cos 


[llx  ' 

-JOjZ 

oo 

lTnx\ 

l~l 

e 

2L  A  cos 
n--0  n 

M 

ja  z 
Jn 


(27) 


where 


°n  * 


k  ;  ?£ 
a 


In  region  B, 


„  /  Hq(x-a)  J^q  ^  Ilf r  (x-a) 

B  cos  (  — -  e  H  r  £_  B  cos  | - - -  e 

I  b  n  0  r  b 


(28) 


whe  re 


'k2  - 
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In  region  C} 


4> 


C 


'wrx 

jv  z  oo 

f^nx  1 

cos  I 

c 

e  f  2_  C  cos 

"Or 

(29) 


where 


Eh 

c 


k  > 


Tin 


jir  n 

2 

,  2 

-  k 

Tin 

—  >  k 


A>  B,  and  C  are  the  amplitudes  of  the  modes  incident  in  regions  A,  B,  and  Cy 

respectively .  The  integers  i ,  qy  and  r  are  arbitrary  „ 

The  coefficients  A  5  B  ,  and  C  are  related  through  the  requirement  that 

n'  n 

the  transverse  E-field  and  H-field  must  be  continuous  across  the  plane  of  the 
discontinuity  (z  ~  0) ■,  Matching  the  transverse  E-field  and  H-field  across 
the  boundary  yields  four  sets  of  equations; 


A  cos 


/£ix 

OO 

*  2„  A  cos 

7i  nx 

— 

B  cosK^0 

1  a  i 

n  0  n 

aJ 

1  b  / 

r  2  B  cos 
n  0  n 


Hn(x-a) 


(30) 


and 


a. A  cos 
i  I  a 


IjUx 


90 


■  2^  a  A  cos 
n-  0  n  n 


TTnx 

a 


■3  B  cos 

q 


^  ,  E  (3  B 
b  n  0  Kn  n 


cos 


f7Tn(x-a) 


(31) 


valid  in  the  interval  c  <  x  <  a 
and  also 

A  cos 


77f  x ' 

OO 

"  2„  A  cos  | 

Vnx 

^  rx 

00 

4  y  r 

cos 

f”n»l 

a 

n  0  n  | 

a 

•  C  cos 

c 

n  -  0  n 

(32) 


and 


a^A  cos 


7I(X 


OC 

^  -  i^nx 

a  A  cos 
n  0  n  n  la 


-  V  C  cos/^4  2  V  C  cos/—' 
r  Ic/nOnn  |c 

I  \ 


( 33 ) 


which  are  valid  in  the  interval  o 
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* 

tor  the  present;  it  shell  be  assumed  that  a/b  is  not  an  integer.  If 
both  sides  of  Equations  (30)  and  (31)  are  multiplied  by  cos  where 

s  is  an  arbitrary  positive  integer;  and  both  sides  of  the  equations  are  then 
integrated  between  the  linits  c  to  a;  two  sets  of  infinite  equations  can  be 


derived .  These  are : 


b(A  $  -  v6:  -  «>■  l 

•  P- 


(l  - 


si°> 


-  <-d*  *-  £ 


*  0=1  a2  -  p3 
n  Ks 


bB  6°  +  £  B  (1-6°)  +  bB6°6°  +  £  B6q  (l  -  6°) 

8  8  2  £  £  '  Q  £  2  J  q  * 


(34) 


and 


- 


*o>6:  +  <-»'  7  ~T~a  <1 


*  P. 


5I>  +  ‘-1>S+1  7  Ji  Vs 


«  V 

n  Ks 


*  bp  B  6°  +  |  p  B  (1  -  6°)  -  bP  b6°6°  -  |  P  B6q  (1  -  6°) 

^s  s  s  2  as  s  sqs  2  s  s  q 


fasi 


(S  ^  0;  1;  2;  0  ->  ®  ) 

By  definition.  A  *  ■=  nA  sin  and  A^  iAain  TLLS. 

*  n  n  a  a 

cb 

The  Kronecker  Delta  b  is  defined  as 

a 


6b  =  1  if  a  -  b 

a 

=  0  if  a  /  b 


A  aim* lar  result  can  be  derived 

c(A6®  ♦  A  )  6°  +  (-1)#+1 

l  os 

=  cC  6  +  —  C 

s  s  2  i 


from  Equations  (32)  and  (33).  Thus^ 


a  2 
a*  -  v 

f  Ts 


/  A ' 

A  ,,  *o„  ,  ,^s+l  V  $  n 

— 1  (1  '  V  "  7nllX“ 


a-  -  \2 

n  8 


(1  -  6°)  +  cC6°6°  +  |  C  6r  (l  -  6°)  (36) 

i  s  r  s  2  s  r 


*  If  a/b  is  a  rational  number,  some  of  the  terms  in  Equations  (34)  and  (35) 
will  become  indeterminate.  The  equations  are  correct;  however;  if  the  in¬ 
determinate  forms  are  replaced  by  their  limits  as  a/b  approaches  a  rational 
value . 
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and 


a  A 

„  cO  ,  vJ-0  ,  ,  .Sir  l  7T  1 

cY  (AO.  -  A  )6  r  (-1)  - — 

o  1  o  s  a  ^2 


a;  -  V2 

I  S 


(1  - 


6;> 


-  _  a  A 

<-i»s  -  f  " " 


a  n  1  a2  .  v2 
n  s 


=cY  C  6°  f  |  V  C  (1  -  6°)  -  cY  C6°6°  -  ~  Y  c6r  (l  -  6°)  (37) 
sss  2  s  s  s  srs  2ss  r 

(s  -=  0,  1,  2 j  . . .) 

Now  refer  to  Equations  (34)  and  (35) .  If  each  line  of  set  Equation  (34) 

is  multiplied  by  P  and  then  sets  (34)  and  (35)  are  added  and  subtracted, 

b 

two  alternate  sets  of  infinite  equations  are  derived .  Thus,  these  can  be 
written  as 


./ 


.  a  .  cOcO  .  ,  „S+1  7T  A' 

bP  Ao  0  +-  (“1)  —  -r - w- 

o  s  f  a  ^  -  p 


(1  - 


_  oo  A 
,  ,,5  ^  ^  n 

('1>  7  nS!  a  .™TT' 

n  s 


|  (3  B  (1  -  6°)  bP  B  6° 
2  s  s  s  o  o  s 


(38) 


and 


oo  A 


.  a  .  cO  ,  .  ,Sfl  7T  A  ,,  £o4  #  ,VS  7T  ^  n 

-bP  A  o  f  (-1)  —  — — —  (l  -  0  )  t  (-1)  ~  £.  ^ — “  IT 

o  o  s  a  a  p  i  a  n-1  a  -  P 

Is  n 


£  P  B  6q  (l 
2  s  s 


6°)  -bP  B6°6° 
q  o  q  s 


(39) 


(S  :  0y  1,  2;  *.) 

Similarly,  from  the  sets  of  Equations  (36)  and  (37),  another  two  sets  of 
equations  can  be  derived  They  are 


CV  A6°6°  - 
o  s  i 


,  _  oo  A' 

<-l>*  1  -  Z 

a  n  1  a 


.  V 
n  s 


|  Y  C  (1  -  6°)  cY  C  6° 
2  s  s  s  o  o  s 


(40) 
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—  »/  ,  _  oo  A' 

,8  "  A'  tov  ,  .  vs+l  V  _  i 


-cVoAo6°  *  .  6°,  ♦  <“l)*+l  ;  &  o-T 

I  s  n 


n  *  \ 


f  Y  C6r  <x  -  6°)  -  cY  c6°6° 

2  s  s  r  *  r  • 


(41) 


(8  s  0)  1>  1)  •••) 


3.2  The  Solution  of  the  Syteii  of  Equations 

In  this  section;  the  scattering  coefficients  are  derived.  This  involves 
the  solution  of  the  systems  of  equations  derived  in  Section  3.1.  Two  separate 
cases  are  considered.  They  are  a)  A  =  C  =  0  but  B  *  0;  and  b)  B  =  C  -  0  but 
A  4  0.  Recall  that  A,  B,  and  C  are  the  mode  amplitudes  of  the  arbitrary  modes 
incident  in  the  three  regions  as  defined  in  Section  3.1.  The  solution  of  the 

systems  of  Equations  (35)  and  (37)  with  A  -  C  =  0  gives  the  elements  of  the 

BB  CB  AB 

scattering  matrices  S  ,  S  ,  and  S  in  terms  of  the  mode  amplitudes  B^, 

BB  CB  AB 

and  A  .  Explicitly,  these  are:  S  ~  -  B  /B,  S  =  -  C  /B.  and  S  =  A  /B. 
m  mn  m  mn  m  mn  m 

The  electric  field  or  voltage  sign  convention  is  followed  when  defining  the 
scattering  coefficients  as  explained  in  the  previous  section.  Similarly,  the 

solution  of  the  systems  of  Equations  (35)  and  (37),  with  B  :=  C  -  0  gives  the 

AA  BA  CA 

elements  of  S  ,  S  ,  and  S  in  terms  of  the  mode  amplitudes  of  that  problem. 

AA  BA  CA 

These  are  expressible  as:  S  -  -  A  /A,  S  =  B  /A.  and  S  =  C  /A. 

mn  m  mn  m  m 

It  is  not  necessary  to  go  through  the  formal  solution  of  a  system  of 

CC  BC  AC 

equations  to  find  the  elements  of  S  ,  S  ,  and  S  since  they  follow  by  a 

BB  CB  AB 

simple  transformation  from  the  elements  of  S  ,  S  ,  and  S  ,  respectively. 

>BB  CB 

If  the  dimensions  b  and  c  are  interchanged  in  the  expressions  for  S  and  S  , 

mn  tan 

CC  BC 

then  the  resulting  expressions  are  identical  to  S  and  S  ,  respectively. 

mn  mn 

Similarly,  if  the  dimensions  b  and  c  are  interchanged  in  the  expression  for 
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S**  and  then  the  total  expression  is  multiplied  by  a  factor  of  (-l)”1.  the  re¬ 
nin 


suiting  expression  is  identical  with  S 


.AC 

mn 


BB  CB  AB 

3»2ol  Derivation  of  the  Elements  of  S  »  S  ,  and  S 

With  A  and  C  set  equal  to  zero,  Equations  (35)  and  ^37)  can  be  written 


as 


oo  a/  abP 

A  s-hr  ♦  <->s  w  ^  '  -nr 


(A  -b6°)6° 
o  q  s 


(42) 


and 


00 

A 


n 


a  _ 
n 


acY 

O  rO 

-z—-  A  6  -r  0 
7T  os 


(s  -S  0,  1.  2 


(43) 


where  q  can  be  any  integer. 

The  above  set  of  equations  can  be  solved  by  the  function-theoretic 

technique.  A  general  discussion  of  the  function  theoretic-technique  is  given 
12 

by  Collin  .  A  meromorphic  function  f(w)  is  constructed  in  such  a  manner 
that  it  will  generate  an  infinite  set  of  equations  which  is  formally  identical 
with  the  original  set  of  equations  when  it  is  integrated  around  the  correct 
contour.  The  form  of  f(w)  depends  on  whether  q  is  zero  or  non-zero  To 
simplify  the  discussion,  consider  first  that  q  -  0,  i.e..,  a  TEM  mode  is  used 
to  excite  region  B. 

A  f (w)  is  desired  such  that 


lim 


*  00 


277  j 


i 


f  (W  ) 

“w"-  T 


n 


n 


^  r<a  > 
oc  n 

S  a  TT 
n  r  l  n  s 


f <P  )6 

o  s 


(44) 


where  r(an)  Is  the  residue  of  f(w)  at  the  pole  w  an^  and 


abP 

f(P  )  *  — 

o  « 


(A  -  B) 

o 


TJ 


(45) 
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and  also. 


lim 

L — *«o 
n 


d^ 

u)  -  V 

s 


=  £ 


n=l 


r(°n) 


+  fcv„>«: 

o  s 


where 


f  (Y_> 
o 


acY 

o 

TT 


A 

o 


(46) 


(47) 


The  contour  Lq  Is  Illustrated  In  Figure  8.  The  function  f(<*>)  is  a  function 
with  simple  poles  located  at 

w  »  o n  «  1,  2,  ...  (48) 

n 

and  simple  zeroes  at 

w  "  Yn,  n  =  1,  2,  ...  (49) 

w  *  Pn,  n  »  1,  2,  ...  (50) 

Furthermore,  by  comparing  Equations  (44)  and  (46)  with  Equations  (42)  and  (43), 
it  is  observed  that 

,  7T nc 

A  a  r (a)/n  sin  — •  (51) 

n  n  a 

The  function  f(<*>)  is  now  determined  to  within  some  integral  function 
p(<*>).  The  function  p(<*>)  can  be  determined  by  examining  the  asymptotic  be¬ 
haviour  of  f(<*>)  as  w— ^oo,  it  can  be  shown  that  in  order  for  the  edge  condi- 
4 

tion  to  be  satisfied. 


f(w) 


o 


-1/2, 


(52) 


excluding  the  poles  and  zeroes  on  the  negative  imaginary  axis. 


Figura  8* 


Location  of  polos  of  f(<*>)  and  the  contour  Lfi  in  tha 
coaplex  w-plans. 
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From  the  above,  one  result  Is  immediate .  Noting  that  P  -  V  =  1c, 

•  o  o 

Equations  (45)  and  (47)  yield  Aq/B  -  b/a.  By  definition,  one  can  write 

S*®  -  b/a  (53) 

oo 


The  function  f(«)  can  be  constructed  in  the  following  manner.  A  function 
f(w)  with  the  necessary  poles  and  zeroes  can  be  written  as 


,  v  ,  n(w  P)  Il(w  V) 
«">  ■  "<“>  ft(<L,a)  1 


(54) 


where,  for  instance. 


n(«,P)  -  Jk  cf»n  -  «>  /jg) 


iSj; 


7  23  4 

The  notation  is  that  used  by  Hurd  *  and  Whitehead  .  The  inclusion  of  the 
exponential  terms  in  the  infinite  product  assures  the  uniform  convergence  of 
the  individual  products.  The  function  p(w),  of  course,  is  the  integral 
function  mentioned  above . 


TTn 


ffn 


Observe  that  asymptotically  for  large  n,  P  ~  -j-r-,  a_  -j-*— ,  and 


n 


n 


Wn 


a  ~  It  can  be  shown  that  if  K(<*>)  is  a  slowly  varying  function  of  w 

n  a 


then  one  can  write 


Vic  tob 

n  ,  wc  *  u>b  Jn7f 

-’(1  "  j  Vn)e  -A(1  "  J 


f(w)  =  K(w)p(w) 


n=l 


Vn‘ 


qo 


«5i 


<*>a 

w®  J  xflh 


(55) 


The  function  f (w)  can  be  written  in  terms  of  the  Gamma  function.  First  recall 
24 

the  identity 


nn  o  +  £> 


u 

n 


-Vu 


(56) 


where  u  is  a  complex  variable,  excluding  the  negative  integers,  V  is  "Euler's 
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Constant’',  and  T(u)  is  the  Gamma  function  of  argument  u  With  the  aid  of 
this  identity,  Equation  (55)  can  be  written  as 


f  (w) 


p(w)K(u>)  77 a 

J  Z  wh  _  ”uc" 

w  b  c  rc-j-jp) 


(57) 


The  function  f(u>)  is  now  written  in  a  form  convenient  for  the  examination 

of  its  asymptotic  behaviour  for  large  values  of  w 

24 

Stirling's  formula  for  the  asymptotic  behaviour  of  the  Gamma  function 
for  large  argument  is  given  by 


r(u) 


(277  ) 


1/2  (u  -  l/2)in  u  -  u 


u- 


>00 


(58) 


This  is  valid  everywhere  in  the  complex  u-plane  except  in  the  vicinity  of  the 
negative  real  axis.  Using  Stirling's  formula  the  asymptotic  behaviour  of 
f (w)  can  be  found  to  be 


f  (w) 


ja  1/2  p(^»K 
2bcj  ^1/2 


(59) 


excluding  the  poles  and  zeros  of  K  Thus,  in  order  for  f(w)  to  exhibit  alge 
braic  growth  at  infinity 


P 


(oj)  oc  e 


Wa 
~jtT 


(60) 


The  constant  of  proportionality  is  determined  by  setting  ftP^)  -  “P  B 
The  relationship  follows  from  Equations  (50)  and  (51)  Finally,  f  (<*<)  can  be 


written 


(|3  _u> 
o 

77 


f  fu} 


bcP  n<u-  6)  n<w  \)  H(p  Q>  -j 

o  „  '  o' 

77  B  ri7F7F>  IT(P  >)  'ii(w>)  e 

o'  o' 


(61) 
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C  and  B  can  be  found  from  f  (w)  .  Note  that  f(w)  satisfies 

9  3 


oo  r(a  ) 


iim 

L - 

n 


2*  j 


f  -  A  r-V  ♦  «-P  ) 


~  n=l  a  +  P 
n  ns 


(62) 


n 


and 


Urn 

L - >oo 

n 


2»Tj 


/ 


»“>  *>  .  2  r<a,|)  +  f(-v  , 

W  +  Y  "  n=l  *  +  V  *'  V 
s  ns 


(63) 


n 


Lr  is  the  same  contour  as  befor  .  Equations  (62)  and  (63)  are  foiptally  iden¬ 
tical  with  Equations  (37)  and  (39).  It  follows  that 


B.  -  <-»*  rh  «-P.> 

s 


(64) 


and 


c  =  (-l)’"1  X(-Y  ) 

s  t  ac  s 


(65) 


where 


2  if  s  :=  0 


1  if  s  ^  0 


It  can  be  shown  that  for  other  modes  of  excitation.  Equations  (64)  and  (63) 


are  valid  providing,  of  course,  the  proper  f (w)  is  used  in  the  calculations. 

BB  CB  AB 

S  ,  S  ,  and  S  follow  immediately  from  Equations  (61),  (64),  and  (65) 
po'  po'  po  * 

and  the  definitions  of  the  scattering  coefficients.  Thus,  for  all  p, 

8B  <-i)p2Pr ,  n<-p  p,  n(-p  v)  n(pc,a,  -j  aL 


90  «  P„  .  n<p  ,p>  n<p  ,y)  n<-p  ,a> 

p  p  o'  o'  p' 


(66) 


and 


S 


BC 

po 


<-l)prl2p  b  I1<-y  .P>  n<-Y  v)  iI(P  a)  -j 


(P  ♦Y_) 


€  Y  a  «UP  P)  n(p  Y>  lii-V  {2) 
p  p  o  o  p 


e 


(67) 
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and  for  p  greater  than  zero. 


.BA 

» 

po 


-bop  Il(a  p)  n ( a.  y)  FI(P  a> 

o  p'  p  o' 

*p  sin^f^  n<P0,P>  n<P0,v> 


(|3  -a  ) 

O  p 

-J  — n —  ai 


(68) 


Q) 


Here  the  superscript  p  in  the  infinite  product  means  delete  the  term  (a  -a  ), 

(  ia  k  /„n  n  p 

n  p,  from  the  product.  By  definition,  a  a 

la 


Consider  now  the  case  of  a  higher  order  TM 


in  region  B  used  to  ex¬ 


cite  the  bifurcated  waveguide.  The  procedure  for  the  solution  of  Equations 
(i2)  and  (43)  is  essentially  the  same  as  before,  differing  only  in  its  details 
First  another  meromorphic  function  of  the  complex  variable  w  is  constructed. 


Call  this  function  f^(w)„  This  is  a  function  with  simple  poles  located  at 


w  -  a  n  1  2  ...  (69) 

xi‘  1  * 

and  simple  zeros  at 

«  - 1,  2,  ...  (70) 

and 

“  -  Pn>  n  »,  1,  2,  ...  (71) 

except  at  u  P^  »  Furthermore, 

f(P  )  (-IT  §£  P  B  (72) 

q  q 

Again,  the  asymptotic  behaviour  of  the  function  at  infinity  must  be  studied 
and  in  order  to  satisfy  the  edge  condition,  f^(w)  ^^(u>  *  as  w  — — ^  oo . 

vote  tha*  in  this  case  has  a  zero  at  ■  p  e  This  is  because  a 

'  1  o 

TEM  mode  is  *ot  excited  in  region  A  by  a  higher  order  TM  excited  in  region  B , 

i  ^p.  A  O  Pus  follows  immediately  if  the  concept  of  reciprocity  is 
'  o 

applied  If  a  TEM  mode  is  excited  in  region  A  no  higher  order  modes  will  be 
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scattered  or  reflected  by  the  bifurcation.  This  is  because  the  incident  TEH 
•ode  is  already  a  normal  node  of  the  syaten.  Hence,  if  no  higher  order  node 
is  excited  in  region  B  by  an  incident  TEH  node  in  region  A,  by  reciprocity, 
no  TEH  node  vill  be  excited  in  region  A  by  an  incident  higher  order  node  in 
region  B.  By  definition,  for  q  >  0, 


3*®  = 
oq 


0 


(73) 


The  proper  f  .(<*>)  can  be  shown  to  be 

(P  -w) 

<-l)qba  P>  P  <w-P  )  n(w  P)  n<u  Y)  n<p  ,a>  —3 -  a 

*,<«)= - a — V, - a —  e  *  L  (74) 

1  aff(0  -w)(P  -p  )  niq;(P  ,p>  n<p  ,y)  n(u>,a) 

q  q  o  q7  q' 


Prom  Equations  (64),  (65),  and  (74)  and  the  definitions  of  the  scattering  co¬ 
efficients,  one  finds  for  q  >  0, 
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m 


.  (P  +P  ) 

(-i)p4,q+1Pq(Pp^Po)  ii(-pp,p)  n<-Pp,Y)  n<pq,q)  -j 

€„P„(P +P  )<P  -P  )  n<q)<p  p)  n<p  y)  n<-p  ,a) 

pppqqo  q'  q'  p 
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CB 
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pq 
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(-Dp4qbPq(Yp^P0)  n<-Yp,P)  n(-Yp,v)  n(pq,q)  »L 

n<q)<P  n<Po ,y)  n*-Y  a) 

pppqqo  q  q  q 


and  for  p  >  0,  q  >  '1, 
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AB 

J 

pq 


<-l)qcaP  <a -P  >  n<a P)  n<a  Y)  n(P *) 
*  P  °  P  P _ 3- 


z*p  sit/^C  (P  -a  xP  -P  >n(q)<p  ,P)n<p  Y>n(p>(tt  ,o) 
\  •/  q  p  q  o  q'  q  p 


(77) 


Mote  that  at  the  beginning  of  this  section,  it  was  assumed  that  the  mode 
incident  in  any  of  the  three  regions  was  a  propagating  mode.  The  derivation 
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of  the  expressions  for  the  scattering  coefficients  was  based  on  this  assump¬ 
tion..  It  is  necessary,,  also,  that  the  scattering  coefficients  be  defined  if 
the  incident  mode  is  evanescent.  To  do  this,  it  is  only  necessary  in  the  ex¬ 
pressions  for  the  scattering  coefficients  to  replace  the  appropriate  real 
propagation  constants  representing  the  propagating  modes  by  the  corresponding 
imaginary  propagation  constants  representing  evanescent  modes.  The  expressions 
derived  above  are  tabulated  in  Table  1. 

AA  BA  CA 

3.2,2  Derivation  of  the  Elements  of  S  „  S  ,  and  S 

With  B  and  C  set  equal  to  zero,  Equations  (39)  and  (41)  can  be  written 


as 
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and 
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(S 


,) 


As  already  noted  in  Se  ion  3,2.1,  a  TEM  mode  excited  in  region  A  will 
not  be  reflected  by  the  bifurcation  nor  will  any  higher  order  modes  be  ex¬ 
cited.  Hence,  for  all  p,  one  can  write 


SM  0  (80) 

po 


and  for  q  >  0, 


Furthermore, 


oq 


(81) 


sBA 

OO 


1 


(82) 


TABLE  l 


BB  CB  AB 

ELEMENTS  OF  S  ,  8  ,  AND  S 


s68  ,  [s88  sB,M 

L  po  p<jJ 

(PQ  +  0_> 

_BB  <“l)  2poc  n<-pp^>  n<-Pp^>  noo,«)  -j - r-E- 

Spo  *  €p  fJp  a  IT?^7P>  n(Po,Y)  Il(-Pp,a)  e 


;  <«»  •  5  -■  (!)} 


gBB  _  <-n|m|<-1  Pq(Pp  .  Po)  n<-pp,p>  n<.pB,v)  n(PQ,°)  1^1 1 

Sp<l  "«PVpP  ♦  V(pq '  po)  n<<1><pq’p>  n<VY>  "(-<Va> 6 


H  *  *  *»  (s)} 


1,0  Vo» 

p  p 


SCB  m  !>  SCB] 

L  PO  pqj 

-»P>‘  aPob  n(-Yp.P)  n<-Yp,v)  no^.a)  .j (i)  ♦  £  »„ 


n<P0,P>  n<Po,Y>  n<-Yp,«> 


(Cb  <-uptq  b  Pq(Yp  P„)  n(-v.,P)  n(-Yp,Y)  n<pa,») 

”  "  S  c  V\  *  P,)<Pq  •  Po>  n<,><P,’P>  "‘V*  n<->p.<‘> 


TT^f  (s)  *  5  (if 


-  H  •  i  -  ft)} 


i 
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and 


iCA 

oo 


(83) 


Only  if  a  higher  order  TM  mode  in  region  A  excites  the  waveguide,  does  the 
system  of  Equations  (78)  and  (79)  have  non-trivial  solutions. 

Again,  the  function-theoretic  technique  is  called  upon  to  solve  the 
above  equations.  A  function,  say  f0(w),  is  constructed  such  that  when  it  is 
integrated  along  contour  it  generates  a  system  of  equations  which  are 
formally  identical  with  Equations  (78)  and  (79),  I/O,  fg(w)  is  *ound  to 


,  (w+a  > 

-A  (w-PQ)  ri(u,p)  n(w^y)  Il(-a^a)  j  — -i-aL 

f2(W)  "  (a  f(3  )(w+a  )  |l(-a  fl)  rt(-Q  ,Y)  rt(«,a)  ® 

I  of  V  * 


(84) 


A  t  is  the  residue  of  f_(w)  at  ^  -  °  so  that 
P  2  p 


pq 


/q  Sin  ia£  \  (Q  .(3  )  rt(a  ,p)  n(a  ,Y)  n(-a  a) 
i  ’  »  I  pp  P  P'  q" 


a 


P  sin 


TTpc 


(a  .p  )  (a  .a  >ri(-a  p)Il(-a  Y)n(p)(a  a> 
q  o  p  q  q  q  p 


(a  ra  ) 
p  q 

J'  ~  "~7T  'aL 

e  (85) 


for  p.  q  >  0* 

Furthermore,  using  Equations  (64)  and  (65),  one  finds  for  p,  q  >  0, 

,  _  <Q  -p  > 

BA  «.i)p+1»q  sin  fPB+p0)n<'pD'p)n<'Pp’v)n('a-r“)  J"~ iifJ>-aL 

Spq  ab(a  *F)'('«'-P  ,3)n(-a--,5>rf(-.p';S) - i -  e  <86) 

q  o  q  p  q  q'  p; 
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.AC 

> 

pq 


(-l)p2/Tq  sinf^4  (V  *P  )II(-V  P)fl(-Y  ,v)ll(-a  ,a> 
_  \  a  /  p  o  p' _ p; _ q 

ac(a  +P  )(a  -Y  )fl(-a  'fP)H(-a  ,Y)II(-v  a) 
q  o  q  p  q  q"  p 


(a  _v  ) 

q  P 

J - tr -  at 


(87) 


A  4  BA  CA 

The  elements  of  S  '  ;  S  and  S  are  tabulated  in  Table  2  As  before. 


to  get  the  scattering  coefficients  for  evanescent  incident  modes,  simply  re- 


TABLE  3 


AA  RA  C.A 

BLEMBtrS  OF  3  ,  S  ,  AND  S 


TABLE  2 

ELEMENTS  OF  8**,  S84,  AND  SCA 


>u 

© 
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place  the  appropriate  imaginary  propagation  constants  by  the  corresponding 
real  propagation  constants  in  the  given  expressions 
3.3  A  Vote  on  the  Numerical  Compute* i  ions 

The  exact  expressions  for  the  various  scattering  coefficients  are  listed 
in  Tables  1  and  2 .  However,  as  the  expressions  stand;  they  are  not  in  a  form 
convenient  for  purposes  of  calculation;  even  using  a  digital  computer.  An 

infinite  product  such  as  n  (P  f  p  1  e  **  where  P  -  -J 

n- 1  n  p  **n  f 

is  very  slowly  convergent ,  especially  for  large  p.  However,  it  is  possible 

to  express  the  various  scattering  coefficients  in  a  form  suitable  for  calcu- 

lations^  using  a  digital  computer. 

BB 

As  an  example,  consider  (see  Equation  (75)).  Assume  that  the  dimen¬ 

sions  of  the  guide  are  such  as  to  allow  only  the  dominant  mode  to  propagate 

in  each  of  the  three  regions  A,  B,  and  CA  It  was  indicated  in  Section  3.1 

BB 

that  it  is  possible  to  write  S  m  terms  of  the  Gamma  Function.  To  explain 

pq 

more  fully,  consider  the  product  above,  which  is  a  term  in  the  expression  for 
Bf3 

S'.  It  is  possible  to  write 
pq 


00 

A 


(P  fp  >  iP  e 

n  p  7Tr* 


bP 

,  n 


11-1 


^  J 


P  b 

p  b) 

l-Hk 

1  7T  n 

r 

(88) 


Using  Equation  (56'  in  Equation  (88),  one  gets 


OO 


n 

n 


1 


1 


OO 

n 

n  1 


(89) 
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It  ci»'  be  shown  that  the  infinite  product  on  the  right  side  of  Equation 
(89)  is  convergent.  Furthermore,  it  converges  much  more  rapidly  than  the 
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original  product*  The  Gamma  Function  ci"  be  calculated  on  the  digital  com¬ 
puter  probably  using  an  already  available  library  subroutine 


It  follows  that  S  can  be  written  as 

pq 


(90) 


00/ 

„5l 


!  .  Jjlli 

1  lift 

IETVEF 


Here  the  prime  means  omit  the  factor 
infinite  product. 


n  q;  from  the 
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4  "HE  ’ VHOMOGFVt.O'S  E-PlANb  BIFURCATION  IN  A  PARALLEL  PLATE  WAVE GUIDE. 

*  -  tV-i-  chap’f  ry  the  generalized  scattering  matr  technique  is  used  to 
derive  *he  solution  to  the  first  of  the1  three  boundary  value  problems  dis¬ 
cussed  in  this  thesis,  viz  ,  the  inhomogeneous  E-plane  bifurcation  in  a 
parallel  plate  waveguide  Refer  to  Figure  2  The  auxiliary  problem  is  the 
boundary  value  problem  associated  with  the  bifurcation  in  a  parallel  plate 
waveguide.  The  bifurcated  waveguide  is  modified  by  introducing  a  dielectric 
slab  with  a  relative  dielectric  constant  in  region  B  The  slab  completely 
fills  region  B. 

With  reference  to  Figure  2,  let  a  TEM  mode  be  incident  from  region  A 

The  scattered  and  reflected  fields  can  be  represented  in  terms  of  TM  modes. 

no 

As  in  the  previous  section,  the  non-zero  field  components  car.  be  derived  from 
a  scalar  function  identical  to  H  Call  this  function  ^(x  The  non-zero 

y 

field  components  can  be  written  as 
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In  regio”  B,  4^  satisfies 


4j  also  satisfies  the  boundary  conditions 


ddi 

0>  x  0,  a,  all  z  and  x  -  c  z  >  0 
and  at  the  dielectric-air  interface^ 

4,(0  )  ;  ^(O*),,  z  ~  0,  c  <  x  <  a 


(96) 


(97) 


(98) 


34<(o  ) 

5z 


94>(o  )  _ 

i nr  »  z.o,  Cl_x<a 


(99) 


also  satisfies  the  edge  condition  at  the  edge  of  the  bifurcation,  i  .e., 

|v4»|  O  d - >  0  (100) 

where  d  ^  [z^  ■?  (x  -  c.)^ 

One  recognizes  that  the  edge  condition  stated  above  is  the  same  as  the  edge 

condition  for  a  bifurcation  without  the  dielectric  present  in  region  B. 

Figure  9  illustrates  the  general  case  of  a  metallic  wedge  with  a  dielectric 

wedge  situated  next  to  it  »  Using  the  condition  that  the  electromagnetic 

26 

energy  density  must  be  integrable  over  any  finite  domain^  Meixner*  shows  that 
an  admissible  singularity  in  |v4*|  at  the  edge  of  the  composite  wedge  is 

given  by  1^4*1  0(d  as  d - *  0,  where  *1  --  ^ .  If  £  0>  "H  -  1/2. 

It  is  possible  to  show  that  each  individual  term  in  the  multiple  scat¬ 
tering  process  as  expressed  by  Equation  (1)  of  Section  2d  satisfies  the 

th 

above  stated  edge  condition,  i.e.,  the  n  par  dal  wave  transmitted  to 
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region  A  satisfies  the  edge  condition  giver  by  Equation  (100)  The  mode  co¬ 
efficients  cf  the  eigenfunction  expansion  ol  the  n**1  partial  wave  in  region  A 

AB  -r-B  BB.  Bn  BA—  — 

are  giver  by  the  elements  of  the  vector  S  1  fS  i  )  S  a  The  vector 
associated  with  the  total  field  reflected  in  region  A  is  the  sum  of  the  in¬ 
finite  number  of  multiple  reflections  That  the  order  of  the  singularity  at 
the  edge  of  the  bifurcation  remains  unaltered  by  the  addition  of  the  infinite 
number  of  terms  in  the  series  given  by  Equation  (1)  has  net  beer  formally 
proved.  However,  the  convergence  of  the  series  given  by  (D  was  proved  in 
Section  2.2. 

It  should  be  noted  that  the  usefulness  of  a  .solution  to  a  problem  is  not 

necessarily  dependent  on  whether  or  not  the  edge  condition  or  any  of  the 

2? 

boundary  conditions  are  exactly  satisfied  Schelkunoff  has  noted  that 
nearly  correct  calculations  are  sometimes  possible  from  solutions  which  only 
to  a  crude  approximation  satisfy  the  boundary  condition:  of  the  problem 

Without  any  loss  of  generality,  it  can  be  assumed  that  the  amplitude  of 
the  incident  TEM  mode  is  unity.  Then  the  total  field  4>  in  region  A,  denoted 
by  4^  can  be  written  as 


Jaz  oo 


— jn  2 

4»  e  °  -  R  e  °  f  2  A  cos  l — -|e 

A  A  n~l  n  \  a 


Wnxl  JV 


(101) 


and,  denoting  the  total  field  in  region  C  by  ^ 


4>  t  e 
C  C  A 


JV  •  r  (lTnx\ 

f  2  C  cos  1 - 1 

n  l  n  \  c  I 


-  i\  Z 

“  n 


(102) 


The  propagation  constants  an  and  are  defined  in  Section  3.1.  R^  is  the 

voltage  reflection  coefficient  for  the  dominant  mode  and  T  .  is  the  trans- 

CA 

mission  coefficient  from  region  A  to  region  C.  The  total  field  in  region  B, 
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denoted  by  /  can  Lc  written  as 
B 


4»  '  Tn» 

B  BA 


oo 


2.  B  cos 
n  1  n 


rn (x-a ) 


-  lP  z 
n 


(103) 


where 


T  is  the  transmission  coefficient  for  the  TEM  from  region  A  to  region  B, 

On 

The  reflectance  and  transmission  matrices  which  characterize  the  dielec¬ 
tric  load  can  be  derived  by  appealing  to  simple  waveguide  theory.  They  are 
both  diagonal  matrices.  The  diagonal  elegants  of  the  reflectance  matrix  are 
given  by 


B  8  ,  ~  8n 

P"  g/n  f  gn 

where  g  ^is  the  wave  impedance  of  the  TM 

n  no 

guide#  height  b,  filled  with  dielectric  and 

TM  mode  m  the  same  waveguide  filled  with 
no 


(104) 

mode  in  a  parallel  plate  wave- 
g^  is  the  wave  impedance  of  the 
air.  Specifically,, 


p' 

gn  oj€  A 


(105) 


and 


(106) 


Call  this  reflection  matrix  1^.  The  transmission  matrix  is  simply  related  to 
the  reflectance  matrix  Let  the  t r»  .lsmission  matrix  be  denoted  by  "36^. 
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Then 


$ 


B 


-B 


D 


t  1 07  i 


,B 

The  general  diagonal  element  of  is  given  by 


B  2< 


V  g  •  g 

*n  6n 


.108) 


Using  Equation  (2>  it  is  now  possible  to  write  the  node  coefficients  in 

Equations  (101)  in  terms  of  known  quantities „  Let  S^n  be  defined  as  the 

T 

vector  (R  -  A  ^  -A  ,  ,  <»)  „  Then 

A  1  A 


in 


sAB  rD  °  '  S 


BB  rB,-l  SBA 


(109) 


BA  BB  AB 

The  elements  of  the  scattering  matrices  S  ;  S  and  S  are  tabulated  in 
Tables  1  and  2.  Equations  (101)  and  (109)  completely  determine  the  fields  in 
region  A , 

/  T 

Let  S_  be  defined  as  the  vector  (T  .  C.  C  ,  ,,.)  »  Then 

CA  GA'  1'  2"  " 


s^r  -  s06  r®  <i  -  sBV  sp^ 


(110) 


Equations  (102)  and  (110)  describe  the  fields  in  rtgion  C= 

Let  S*  be  defined  as  the  vector  (T  . ,  B,  B  ,  Then 

BA  BA  1'  2 


5^  (I  -  S 

PA  D  V  D 


BBrV‘  sBAI 


am 


Equations  003)  and  (111)  describe  the  fields  in  region  B 

Vote,  however>  that  the  order  of  the  matrix  (!  -  sBBr3  is  infinite  ,  and 
no  method  is  now  known  to  invert  this  matrix  exactly,  Tt  will  be  shown  that 


accurate  calculations  of  such  desirable  quantities  as  the  reflection  co- 
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efficient  R  using  Equation  (109),  arc  possible  by  working  with  finite  order 

matrices,  i  e  ,  by  including  in  the  calculation  a  finite  number  of  elemen  s  in 

BBrB 

the  matrix  (1  -  S  i  >  The  rapid  convergence  of  the-  solution  with  an  increase 

in  the  truncation  size  is  demonstrated  in  the  following  discussion 

BB  P 

Consider  Equation  (1091  Let  the  matrix  (I  -  S  r  )  be  truncated  to  a 
mi'rix  of  finite  order  N(  N  -  1  ;  This  means  that  only  the  first  N  rows  and 
c  lunr-s  of  the  truncated  matrix  (I  -  are  included  in  the  calculation 

Let  the  determinant  of  the  truncated  matrix  be  denoted  by  Let  the  deter¬ 

minant  of  the  minor  of  the  truncated  matrix  obtained  by  striking  out  the  *irst 
row  ar.d  column  be  denoted  by  Aj  ^  ^  ^  .  Then,  it  can  be  shown  that  is 
given  approximately  bv 


b  B  A11  (N- 1 1 


(112) 


Oi  course,  if  V  1,  then  Equation  (109)  reduces  tc  a  scalar  equation  for  R 

ri 

which  car  be  written  as 


b 


B 

Po 


BB  B 
S  p 
oo  r’o 


where 


sDB 

OO 


is  given  by 


Po  as  defined  above  is  giver  by 

B  1  -Jk 


(113) 


'114' 


(llo) 
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It  was  demonstrated  numerically  that  the  approximate  expression  for 
given  by  Equation  (1121  rapidly  converges  to  a  limit  as  the  order  v  of  the 
truncated  matrix  is  increased-  The  results  for  a  particular  example  are  shown 
in  Table  3  below.  The  choice  ot  parameters  for  this  set  of  calculations  are 
K  2.5>  a/\  =  .339,  and  c/a  .5.  Polystyrene  has  a  relative  dielectric  con¬ 
stant  around  2.5  and  is  a  commonly  used  dielectric  The  actual  values  of  the 
elements  comprising  the  first  five  rows  and  columns  of  the  matrix  (I  - 
for  this  set  of  parameters  used  in  this  calculation  are  shown  in  Table  4. 

The  results  of  the  calculations  indicate  tha*  che  major  contribution  to 


the  value  of  the  reflection  coefficient  R  comes  from  the  terra  given  by 

A 

Equation  (1A3),  i  .e  the  term  due  to  the  TEM  mode  alone  =  The  contribution 

to  the  value  of  RM  due  to  the  higher  order  TM  modes  is  small.  Indeed  for 

A  r,o 

a/X.  «  .5,,  Equation  (113)  is  an  accurate  expression  for  R  - 


TABLE  3 

Reflection  Coefficient  R.  for  I nhomogeneous  E-Plane 

A 

Bifurcation  (K  2.5,  a/X  -  .339,  c/%  .5) 


Rank  If  of  the 

Truncated  Matrix 

(i  -  sEBr£ 

Reflection  Coefficient 

*A 

1 

--lOSeJ510 

_  o 

2 

-.107ej54 

-  o 

3 

-,107.J5-4 

o 

4 

-  107ej5c4 

A  system  of  infinite  order  linear  algebraic  equations  can  be  derived  for 

the  inhomogeneous  bifurcation.  The  derivation  of  thi9  system  of  equations, 

16 

giver  by  Mittra  and  Pace  ,  is  essentially  the  same  as  the  derivation  of  the 


TABLE  4 


ELEMENTS  OP  FIRST  FIVE  ROWS  AND  COLUMNS  OF  (I  - 
<K  =  2.5,  a/\  ---  .339,  c/a  =  .5) 

N  =  1  N  =  2  N  -  3  N  - 

.9436  4-  j  .0973  |  .1133  -  j.0668  -.0762  4-  j .0447  .  0626  -  j .0368 

.0226  4-  j  .0386  .  9685  -  J  .0260  .0307  f  j  .0174  -.0289  -  j  .0143 

(I  -  SBV)  =.  -.0077  -  j  .0131  .0157  +  J  .0089  .9833  -  J  .0059  .0167  4-  j  .0049 

.0042  4  J  .0071  -.0097  -  j  .0048  .0111  4-  j  .0032  .9885  -  j  .0026 

-.0027  -  j  .0046  .0068  4-  J  .0031  -.0080  -  j.0021  .0086  4-  j  .0017 


I 


TABLE  4 


IVE  RCWS  AND  COLUMNS  OF  (I  -  SBBI^) 
5>  a/\  -  .339,  c/a  ~  .5) 


N  =  2 


N  -  3 


N  -  4 


-  j ,0668 

-  J  .0260 


+  J ,0089 


-  j .0048 


- .0762  +  j .0447  .0626  -  j .0368 

.0307  +  j .0174  -.0289  -  j .0143 


.9833  -  j .0059 


.0167  +  j .0049 


.0111  +  j  .0032  .9885  -  J  .0026 


- .0538  + 
.0268  + 
-.0161  - 
.0114  + 


.0316 

.0123 

.0042 

.0023 


+  J.0031  -.0080  -  J .0021 


0086  +  j .0017 


.9914  -  j.0015 
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system  of  equations  corresponding  to  the  homogeneous  bifurcation  (see  Section 
3.1).  However,  unlike  the  equations  derived  ir  Section  3.1,  the  set  of  equa¬ 
tions  for  the  inhomogeneous  bifurcation  cannot  be  solved  exactly.  An  iter¬ 
ation  technique  must  be  used  for  its  solution.  A  relationship  has  been  shown 
between  the  iterative  solution  of  the  system  of  equations  pertaining  to  the 
inhomogeneous  bifurcation  and  the  scattering  matrix  formulation  of  the  same 
problem.  The  Iterative  solution  of  the  above-mentioned  system  of  equations 

is  not  nearly  as  convenient  for  purposes  of  calculation  as  say  Equation  (109). 

28 

Cronson  in  the  only  known  paper  that  deals  specifically  with  the  problem 
of  the  Inhomogeneous  E-plane  bifurcation  has  derived  an  equivalent  circuit  for 
the  Junction.  Refer  to  Figure  10.  Cronson  expresses  the  equivalent  junction 
capacitance  in  terms  of  an  infinite  series  of  sine  terms  with  constant  co¬ 
efficients.  The  coefficients  are  the  solutions  to  a  system  of  infinite  order 
linear  algebraic  equ^ions.  Cronson  finds  that  it  is  not  possible  to  solve  his 
system  of  equations  exactly.  A  method  of  approximation  is  employed.  He  re¬ 
sorts  to  solving  a  truncated  set  of  equations.  Specifically,  he  solves  a 
system  of  equations  of  rank  6.  Hie  use  of  a  digital  computer  is  required  to 
carry  out  all  of  the  computations,  including  the  computations  of  the  elements 
of  the  sixth  order  matrix  to  be  Inverted. 

Cronson  checks  the  validity  of  his  approximations  by  comparing  his  results 

for  the  case  A  =  1  with  the  results  given  by  Marcuvitz  in  the  ’Waveguide  Hand- 

17 

book"  for  the  homogeneous  bifurcation.  Working  with  the  normalized  capaci¬ 
tance  Cy(Cy  ~  -  ) ,  Cronson  finds  the  following  percent  errors  in  his 

calculations;  c/a  -  .5,  a/X  •  .5,  2 .2^  c/a  .5,  a/X  .3,  3.46%;  c/a  -  .5, 

s/X  -  .5,  8.52%. 

The  percentage  error  in  computing  R  from  is  less  than  the  percentage 


error  in  C  itself  .  For  Instance,  for  the  set  of  parameters  a/X  -  .339, 


Figure  10.  Equivalent  circuit  for  inhomogeneous 
E-plane  bifurcation. 


Z  Z  and  Z  _  are  the  characteristic  impedances  of  Regions  A 
oA  oB  oL 

B,  and  C,  respectively. 
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c/a  -  .5  and  K  =  2 .5,  Cronson  computed  -=  .31.  From  this  value  for  C^3  one 

1 5  5° 

calculates  R  =  -,108eJ  *  .  This  is  to  be  compared  with  the  results  shown  in 

*  o 

Table  3.  Using  a  matrix  (I  -  SBBI^>  of  fourth  order,  R^  -  -.lOTe'-’5’4  was 
calculated . 


5f 


5.  THE  E-PLANE  METALLIC  STEP  DISCONTINUITY 
Much  has  been  written  in  recent  years  about  the  problem  of  the  step  dis¬ 
continuity  in  a  waveguide.  Refer  to  Figure  1  .  The  inclusion  here  of  a  dis¬ 
cussion  of  the  step  discontinuity  problem  is  justified  on  two  counts.  First, 
from  the  viewpoint  of  studying  the  generalized  scattering  matrix  formulation, 
the  step  discontinuity  is  particularly  interesting  because  it  represents  a 
'worst  possible  case'.  The  configuration  of  the  auxiliary  problem  is  modified 

by  placing  a  perfectly  conducting  wall  in  region  B  flush  with  the  plane  of  the 

__B 

junction,  which  implies  r  *  -I .  It  is  expected  that  the  effect  of  the  higher 
order  modes  is  greater  in  this  case  than  in  the  preceding  problem  of  the  in¬ 
homogeneous  bifurcation.  It  is  desirable  to  show  that  even  in  an  extreme  case, 
important  quantities  such  as  the  reflection  coefficient  R  can  be  computed 

A 

easily  because  of  the  rapid  convergence  of  the  matrix  series  expansion. 

The  second  point  is  that  the  step  discontinuity  problem  warrants  atten¬ 
tion  for  its  own  sake.  It  is  a  classic  problem,  studied  by  a  number  of  authors 

29  17 

using  more  established  methods.  Macfarlane  and  Marcuvitz  have  found  quasi- 
static  solutions  to  the  step  discontinuity  problem.  Their  methods  differ  in 
the  exact  details,  but  essentially  are  the  same.  The  equivalent  susceptance 
of  the  waveguide  junction  is  formulated  in  terms  of  an  integral  equation,  the 
exact  solution  of  which  is  not  possible  in  general .  The  equation  is  solvable 
for  the  case  of  U  (2 v/\)  —  0.  The  static  field  problem  is  solved  by  simpli¬ 
fying  the  original  problem  through  conformal  transformations.  Extensive  re¬ 
sults  are  tabulated  in  the  "Waveguide  Handbook’  .  Only  the  case  of  single  mode 
propagation  is  considered.  In  this  frequency  range  *he  equivalent  static 
method  yields  very  accurate  answers.  However,  the  method  becomes  quite  in¬ 


volved  for  the  case  when  the  dimensions  of  the  guide  are  such  that  several 
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■odea  propagate  i.e.y  a  multi -mode  waveguide.  The  multi-mode  problem  is 
assuming  Increased  importance.  It  has  applications  in  the  field  of  milli¬ 
meter  wave  propagation30  and  the  study  of  VLF  propagation31. 

32 

Schwinger  has  solved  the  above-mentioned  integral  equation  for  the 
equivalent  susceptance  by  means  of  a  variational  technique  .  However^  the 
accuracy  of  the  variational  technique  is  dependent  upon  the  choice  of  the  trial 
function  and,  in  this  sen se,  is  not  a  deterministic  method.  Furthermore,  the 
choice  of  the  trial  function  in  the  variational  technique  is  not  at  all  straight 
forward  for  multi-mode  propagation. 

22 

As  previously  discussed  in  Section  2 .3,  Williams  has  applied  the  Wiener- 
Hopf  technique  to  the  step  discontinuity  problem.  Williams'  method  was  out¬ 
lined  in  Section  2.3  and  the  differences  between  his  approach  to  the  step 
discontinuity  problem  and  the  approach  based  on  the  generalized  scattering 
matrix  technique  were  discussed.  Williams  also  includes  some  numerical  results 
in  his  paper.  For  the  case  of  a/\  <  .5,  his  numerical  results  are  in  close 
agreement  with  Marcuvitz's  results.  Williams  also  discusses  the  situation  when 
two  modes  are  allowed  to  propagate  in  the  larger  channel  of  the  guide  (region  A) 
including  numerical  values  of  the  square  of  the  magnitude  of  the  reflection 
coefficient  R^.  Williams'  results  will  be  used  for  comparison  with  the  results 
reported  in  the  section. 

There  is  still  active  interest  in  the  step  discontinuity  problem  as 

33 

evidenced  by  the  most  recent  paper  on  the  subject  by  Magnus  and  Fox  .  In 
uiA  paper,  the  problem  is  treated  as  an  infinite  set  of  inhomogeneous  linear 
equations.  They  are  solved  formally  by  a  perturbation  technique. 

With  this  brief  introduction,  now  consider  the  problem  of  the  step  dis¬ 
continuity  in  terms  of  the  generalized  scattering  matrix  formulation.  With 
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reference  to  Figure  1,  assume  that  a  TEM  mode  of  unit  amplitude  is  incident  in 
region  A,  traveling  in  the  positive  z  direction.  As  in  Section  4,  the  non¬ 
zero  field  components  H  ,  E  ,  and  E  can  be  derived  from  a  scalar  function 

y  x  z 

4>(x>  z)  ■  In  this  case,  4*  satisfies  the  homogeneous  Helmholtz  equation  to¬ 
gether  with  the  boundary  conditions 

■gj  =  0,  x  i  0,  a,  all  z  and  x  -  c,  z  >  0.  (116) 

and 

.=  0,  z  -  0  c  <  x  <  a  (117) 

i)z  '  —  — 

4*  also  satisfies  the  edge  condition  given  by 

|v  4>j  -  0<d_1/3),  d - (H8) 

where  d  -  [  (x  -  c)^  r 

Let  the  fields  in  region  A  be  expanded  in  the  cosine  series  given  by 
Equation  (101).  The  mode  coefficients  of  the  reflected  field  are  given  by 


s88)'1 


019) 


-  r,B 

where  S,  is  the  column  vector  defined  in  Section  4.  Here,  1  -T  where  I 
in  * 

is  the  identity  matrix. 

Similarly,  let  the  fields  in  region  C  be  expanded  in  the  cosine  series 
given  by  Equation  (101).  Then  the  mode  coefficients  expressed  by  the  column 
vector  S  are  given  by 


(i 


BB  -1 

r  S  ) 


(120) 
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The  fields  in  region  B,  of  course  are  identically  zero. 

That  the  solution  expressed  by  Equations  (119)  ard  (120)  satisfies  the 
edge  condition  given  by  Equation  (118)  has  not  beer  shown.  In  order  to  do  this, 
one  should  exanine  the  asymptotic  behaviour  of  the  higher  order  mode  co¬ 
efficients.  i.e.,  the  higher  order  elements  of  the  vectors  S  or  S  ..  In 
this  section,  only  the  reflection  coefficient  R  is  computed.  However,  the 

n 

BB  BB  BB 

proof  of  the  convergence  of  the  Neumann  series  I  +  S  +  S  S  +  ...  is 
given  in  Section  2.2.  Also,  it  is  demonstrated  in  this  section  that  for 
0  <  a/X.  <  1.0,  the  computed  values  for  R  are  in  close  agreement  wi*h  figures 

i\ 

computed  from  the  expression  given  in  the  'Waveguide  Handbook",  or  alternately, 
with  figures  reported  by  Williams 

Now  in  Equation  (119),  let  the  matrix  (I  i-  S  )  be  truncated  to  a  matrix 
of  order  N,  N  >  1 .  As  in  Section  4,  let  the  determinant  of  the  truncated  matrix 


22 


.BB. 


be  denoted  by  and  let  the  determinant  of  the  minor  of  the  truncated  mat- 

(“) 

rix  obtained  by  striking  out  the  first  row  and  column  be  denoted  by  A 

ll(N-l) 

Then  Ra  is  given  approximately  by 
A 


b  ^lKN-l) 
r  ~ - - - 

A  *  a  A(N) 


(121) 


If  N  1,  then 


oo 


(122) 


where  Soo  is  given  by  Equation  (114)  for  0  a,  \  <  0  5 

Even  in  this  case,  it  was  demonstrated  numerically  that  the  approximate 

Equation  (121)  for  R.  converges  rapidly  to  a  limit  as  the  order  N  of  the 

A 
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,  BB, 

truncated  matrix  (I  +  S  )  increases  .  The  results  of  two  sets  of  calculations 

are  shown  in  Tables  5  and  6.  The  results  of  the  calculations  indicate  that  in 

BB 

the  range  0  <  a/K  <  .5,  a  truncated  matrix  (I  *  S  )  of  order  four  or  five  is 

sufficient  for  purposes  of  accurate  computation. 

For  c/a  =  .5  and  a/\  ~  .339,  R.  -  -.393e^^^*^  was  calculated  from  the  eA- 

A 

pression  for  the  equivalent  susceptance  of  the  junction  given  by  Marcuvitz  in 

it  121  .3° 

the  "Waveguide  Handbook".  Compare  this  with  the  value  R.  -  .392eJ  *  calcu- 

A 

lated  from  Equation  (121)  with  N  -  5.  The  actual  values  of  the  elements  com- 

BB 

prising  the  first  five  rows  and  columns  of  the  matrix  (I  f  S  )  for  this  set 


of  parameters  used  in  the  calculations  are  shown  in  Table  7.  In  the  second 

i20  9° 

example  with  c/a  -  .326  and  a/X,  -  „3,  R  ■=  -  .215eJ  '  was  computed  from 

A 

'20  3° 

Marcuvitz's  figures.  RA  -  .215eJ  *  was  calculated  from  Equatlan  (121) 

with  N  4  . 


TABLE  5 

Reflection  Coefficient  R  for  Step  Discontinuity 
(a/X.  --  .339,  c/a  =  .5) 


Rank  N  of 

Truncated  Matrix 

(1  ♦  sBB> 

Reflection  Coefficient 

ba 

1 

- .379ejl8°8 

2 

-.388ej2°-7° 

3 

-,391ej21-10 

4 

-.392.j21  2° 

5 

-.392ej21-3° 
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TABU.  6 

Reflection  Coefficent  P.  for  Step  Discontinuity 

A 


(a/X  3,  c/a  -•  326  > 


Rank  of  Truncated 

BB 

Matrix  (I  4  S  ) 

Reflection  Coefficient 

ba 

1 

-w19-6° 

2 

-.214.jl9-8° 

3 

-.2156j20-°° 

4 

-.215,J20-3° 

Next,  consider  the  situation  when  the  TEM  and  TM^  modes  are  allowed  to 


propagate  in  region  A,  but  only  the  TEM  mode  propagates  in  region  C .  The  re¬ 
flection  coefficient  RA  was  computed  using  Equation  (121)  and  Equation  (122). 
The  results  of  the  numerical  calculations  for  two  examples  are  listed  in 

Tables  8  and  9.  One  can  compare  these  results  with  those  published  by  Williams 

12 


For  a/X  -0.7  and  c/a  r  0.5,  Williams  computes  R 


'A 


.28.  Using  Equation 


(121)  with  N  -  5,  one  computes  R 


.29.  The  values  of  the  elements  of  the 


.BB. 


first  five  rows  and  columns  of  (I  +  S  )  for  this  set  of  parameters  are  shown 

.  12 


in  Table  10.  far  a/X  --  0.9  and  c/a  -  0.5,  Williams  computes 

|2 


.25.  Again 


using  Equation  (121),  N  -  5,  one  can  compute 


25,  Thus,  there  is  close 


agreement  between  the  results  reported  by  the  author  and  Williams'  results  even 

when  more  than  one  mode  is  allowed  to  propagate  in  region  A  The  dominant 

BB 

contribution  to  R.  comes  fror  Equatio"  (122>  with  S  given  by 

A  oo 


(123) 


TABLE  7 


ELEMENTS  OP  FIRST  FIVE  ROMS  AND  COLUMNS  OF  (I  +  SBB) 

(a/\  =  .339,  c/a  =  .5) 


N  =  1  N  2  N 

1.2508  -  j .4326  |  -.2411  f  j .1416  .1736  -  j .1019 

-.1006  -  j  .1714  1.0667  4-  j .0551  -.0698  -  j .0396 

.0342  -s*  j  .0583  -.0332  -  j  .0188  1.038  i-  j  .0135 

-.0186  -  j  .  03 16  .0206  4-  j  .0102  -.0252  -  j  .0073 

.0121  +  j  .0205  -.0144  -  j  .0066  .  0183  4-  j  .0048 


3 

-.1445  + 
.0667  - 
-.0386  - 
1.0265  + 
- .0198  - 


N  = 

j .0849 
j .0329 
j .0112 
j .0061 
j .0040 


TABLE  7 


DD 

OP  FIRST  FIVE  ROWS  AND  COLUMNS  OF  (I  +  S) 
(a/\  =  .339,  c/a  =  .5) 


=  1 

N 

2 

N  . 

-  3 

N  : 

=  4 

J  -,2411 

f 

j  .1416 

.1736  - 

j .1019 

- .1445 

+ 

j .0849 

.1249 

- 

1 .0667 

+ 

j  .0551 

- .0698 

-  j .0396 

.0667 

- 

j .0329 

- ,0622 

- 

-.0332 

- 

j  ,0188 

1.038  r 

j .0135 

- .0386 

- 

j  .0112 

.0373 

+ 

.0206 

•t* 

j  .0102 

-  .  0252 

-  j ,0073 

1 .0265 

+ 

j .0061 

-.0263 

- 

-.0144 

- 

j  .0066 

.0183 

+  j .0048 

- .0198 

- 

j .0040 

1.0200 

j  ,0733 
j ,0286 
j ,0097 
j ,0053 
j ,0034 
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TABLE  8 

REFLECTION  COEFFICIENT  R4  FOR  STEP  DISCONTINUITY 

A 

TWO  MODES  PROPAGATING  IN  REGION  A 
(a/X  =  0.7,  c/a  =  0.5) 


Rank  N  of  Truncated 

no 

Matrix  <1  +  8  ) 

Reflection  Coefficient 

ra 

1 

-.52^4-6° 

2 

-.52ej2’8 

3 

-.52ej2,2° 

4 

«.52ej2*2 

5 

-W22” 

TABLE  9 

REFLECTION  COEFFICIENT  R.  FOR  STEP  DISCONTINUITY 

A 

TWO  MODES  PROPAGATING  IN  REGION  A 
(a/X  0.9,  c/a  =  0.5) 


Rank  N  of  Truncated 

BB 

Matrix  (I  +  S  ) 

Reflection  Coefficient 

ra 

1 

-.50eJ2-6° 

2 

-.50.J3-7° 

,  O 

3 

-  ,50eJ  * 

,  _o 

4 

-.50eJ1’7 

5 

-.50eJ1-7° 

TABLE  10 


ELEMENTS  OF  FIRST  FIVE  ROWS  AND  COLUMNS  OF  (I  +  S®B) 
TWO  MODES  PROPAGATING  IN  REGION  A 
<aA  =  0.7,  c/a  =  .5) 


N  =  1  N  =  2  N^ 

.9553  -  j .0762  |  .0164  +  j .1499  .  0218  -  j.0906 

-  .2849  +  j  .0312  1.2251  j  .0817  -.1591  -  j.0450 

.0671  +  j .0161  -.0615  -  j .0174  1.0539  +  j .0087 

-.0332  -  j .0112  .0337  +  j .0083  -.0323  -  j .0056 


3  N  = 

-.0236  +  j .0702 

.1359  +  j  .0337 

-.0505  -  j .0087 

1 .0317  +  j .0044 


.0206  +  j .0079 


-.0222  -  J .0051 


.0225  +  j .0034 


.0228  -  j .0028 


TABLE  10 


DD 

rs  OF  FIRST  FIVE  ROWS  AND  COLUMNS  OF  (I  +  S  ) 
TWO  MODES  PROPAGATING  IN  REGION  A 
(a/X.  =  0.7,  c/a  =  .5) 


N  = 

1 

N  = 

:  2 

N  E 

■=  3 

N  = 

=  4 

lI 

.0164 

+ 

j .1499 

.0218 

- 

j .0906 

-.0236 

+ 

j  .0702 

.0228 

- 

j  ,0593 

2 

1.2251 

f 

j  .0817 

-  .1591 

- 

j .0450 

.1359 

+ 

j  .0337 

-.1215 

- 

j  .0280 

1 

-.0615 

- 

j  .0174 

1.0539 

+ 

j  .0087 

-.0505 

- 

j  .0087 

,0476 

j  ,0074 

2 

.0337 

+ 

j  .0083 

-  .0323 

- 

j  .0056 

1 .0317 

4- 

j  .0044 

-  ,0309 

- 

j  .  0038 

9 

-.0222 

- 

j  .0051 

.0225 

+ 

j  .0034 

-.0228 

- 

j  .0028 

1 ,0226 

+ 

j  .0024 

o> 

w 
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6  THE  TRIFURCATED  WAVEGUIDE 

This  section  is  concerned  with  the  boundary  vilue  problem  associated  with 
two  s«mi~inf inite  plates  in  a  parallel  plate  waveguide.  Refer  to  Figure  3. 

Let  a  TEH  mode  be  incident  in  region  C.  Let  4*  (x,  z)  be  a  scalar  function 
such  that  4*  *  H  .  Then,  E  and  E  can  be  derived  from  4*  using  Equations  (93) 

y  x  y 

and  (94)  with  €  *  .  4*  satisfies  the  homogeneous  Helmholtz  equation  together 

with  the  boundary  conditions 


do* 

■5—  «  0,  x  -  0,  a  for  all  z  and  x  -  c,  h  for  z  >  0  (124) 


and  the  edge  conditions 


Ml  =  O  (d1>2~1/2)'  dl,2 - >°  (125) 

v  ^  r  /  ,2  2 ,1/2 

where  dj  =-  [  (x  -  c)  +  z  J 
and  ~  [  (x  -  h)  +  z  J 

4^(x,z)  can  be  expanded  in  region  C  in  the  cosine  series  given  by 


-•>v  ->V 

e  -  R.  e 


00 

A 


cos 


fr) 


JVnz 


(126) 


where  R^  is  the  voltage  reflectance  coefficient  for  the  TEM  mode.  The  ampli¬ 
tude  of  the  incident  TEM  mode  in  region  C  is  assumed  to  be  one. 

An  expression  for  the  coefficients  of  the  Fourier  series  expansion  given 
by  Equation  (126)  can  be  written  in  terms  of  the  scattering  coefficients  of 
the  auxiliary  problem  and  the  load  in  region  B,  Thus,  one  can  show  that 


,  scc  7  ►  sCBrB  <1  -  s^V1  sK  7 


(127) 
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-//  X  —  T  B 

where  ;  (-R  ,  ^1'  ^23  anc*  c  *  (1,  0,  0*  •  ■>•)  *  The  elements  of  I 

in  this  case  are  Known .  As  mentioned  before>  the  bifurcated  waveguide  is  modi¬ 
fied  by  placing  a  second  semi-infinite  plate  in  region  B,  Thus,.  the  derivation 


r8 


_  B 

of  the  elements  of  1  for  this  problem  is  the  same  as  the  one  followed  in 

AA 

solving  ter  the  elements  of  S  in  Section  2.  To  derive  the  elements  of 

AA 

from  the  expressions  for  S  given  in  Table  2  simply  replace  a  by  b,  b  by  h  -  c. 

mn 

and  b  by  a  -  h . 

Note  the  semi-infinite  plate  in  region  B  is  coupled  electro-magnet ically 

to  region  C  by  the  higher  order  n*no  modes.  A  TEII  mode  scattered  in  region  B, 

traveling  in  the  positive  z  direction,  will  not  be  reflected  by  the  septum  in 

that  region.  The  higher  order  TM  will  be  reflected,  however-  and  will  con- 

no 

tribute  some  to  the  final  value  of  the  reflection  coefficient  R  .  The  dominant 

cc 

term  in  the  series  expansion  for  R^  is  given  by  since  as  shown  by  actual 

calculation,  the  contribution  of  the  higher  order  TM  modes  is  numerically 
3  no 

CC 

small  compared  to  S  .  This  is  to  say  that  the  second  plate  can  be  introduced 
in  region  B  without  appreciably  affecting  the  reflection  coefficient  R  . 

The  reflection  coefficient  R  can  be  computed  from  Equation  (127)  if  first 
each  of  the  matrices  appearing  in  the  equation  are  replaced  by  a  matrix  of 
order  N.  It  was  demonstrated  by  means  of  actual  calculations  that  in  the  case 
of  both  the  step  discontinuity  and  the  inhomogeneous  bifurcation,  the  ex¬ 
pression  for  the  reflection  coefficient  rapidly  converges  to  a  limit  as  the 
order  N  of  the  truncated  matrices  increases.  As  examples,  the  numerical  re¬ 
sults  of  two  sets  of  calculations  are  cited  in  Tables  11  and  12,  The  parameters 
used  in  the  first  example  are  h/a  .5,  c/h  .5,  and  a/\  -  .4 .  In  the  second 

example,  they  are  h/a  -  .326 ,  c/h  .326,  and  a/\  ,3.  In  both  of  these 

examples,  the  two  plates  are  asymmetrically  situated  with  respect  to  the  center 


TABLE  11 


REFLECTION  COEFFICIENT  R  FOR  TR’ FURCATED  WAVEGUIDE 

v 

(h/a  -  .5,  c/h  .5/  a/\  *  .4) 


No.  of  Modes  (N) 

Included  in  Calculations 

Reflection  Coefficient 

Rc 

1 

.SOOe-J75-®0 

2 

.470eJ78-4° 

3 

.470e-J78'9° 

4 

.475a-J77-9° 

6 

•474e-J78'10 

TABLE  12 

REFLECTION  COEFFICIENT  R  FOR  TRIFURCATED  WAVEGUIDE 

v 

(h/a  =  .326,  c/h  «  .326,  a/\  =  .3) 


No.  of  Modes  (N) 

Included  in  Calculations 

Reflection  Coefficient 

Hc 

1 

.326e-J46-9° 

2 

•  327e”^48  *9° 

3 

.324e-J47'5° 

4 

.324e‘J47-5° 

6 

•324e-J47-4° 
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line  of  the  waveguide  described  by  x  a/2:  It  is  shown  below  that  a  simple 

expression  can  be  derived  for  R  if  the  two  plates  are  symmetrically  placed 

0 

with  respect  to  x  -  a/2, 

It  is  possible  to  formulate  the  problem  of  the  trifurcated  waveguide  in 

34 

terms  of  a  system  of  Wiener-Hopf  integral  equations  A  E  Heins  discusses 
the  special  case  of  an  arbitrary  number  of  equally  spaced  semi-infinite  plates 
in  a  waveguide.  The  set  of  integral  equations  can  be  formulated  in  terms  of 
the  unknown  current  densities  on  each  of  the  semi-infinite  plates  For  the 
case  of  the  trifurcated  waveguide,  the  system  of  infegral  equations  are  of  the 
form 

00 

2  r  /  /  / 

I  /  K.  (z  -  z  )  J.  (z  >  dz  F  <z>  0  (128) 

J  -  l'/0  U  1  J 

for  z  >  0  and  i  .  la  2,  K.  (z  -  z  \  are  linear  combinations  of  the  Green’s 

*  iJ 

functions  used  in  formulating  the  integral  equations  and  F  f z 1  is  the  form  of 

J 

th 

the  propagating  modes  in  the  j  duct.  The  solution  of  Equatior'  (128)  is  com¬ 
plicated,  lb  solve  the  above  set  of  equations,  the  Wiener-Hopf  technique  must  be 
generalized.  It  is  necessary  to  factorize  the  determinant  of  the  matrix  whose 

elements  are  the  Fourier  transforms  of  the  kernels  K.  (zt,  Ir,  this  instance, 

^  J 

one  must  factorize  a  determinant  of  order  two,  Heins  has  discussed  this  pro¬ 
blem  in  general  terms,  but  he  has  not,  as  far  as  it  is  known,  published  the 

35 

actual  solution  to  the  problem.  Wu  and  Wu  in  a  paper  published  much  later 
than  Heins  state  that  in  the  case  of  coupled  Wiener-Hopf  integral  equations, 
solutions  are  not  known  except  for  the  cases  where  reduction  *o  *  single 
equation  is  possible 

The  exact  solution  to  the  trifurcated  waveguide  is  possible  when  the 
plates  are  symmetrically  spaced  with  respect  to  the  center  line  x  a,  2  as 
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mentioned  above.  This  is  discussed,  next. 

It  is  desired  to  find  the  reflection  coefficient  R^  for  the  dominant  mode 
incident  in  region  C,  Refer  again  to  Figure  3.  This  mode  of  excitation  can 
be  broken  into  two  separate  cases  of  even  and  odd  excitation.  Consider 
Figure  11.  The  area  of  the  three  smaller  ducts  to  the  right  of  the  plane  of 
the  Junction  are  labeled  B',  B"  ,  and  C.  In  the  case  of  even  excitation 
(Figure  11a),  regions  B'  and  C  are  excited  in  the  TEM  mode  so  that  the  incident 
fields  are  in  phase.  In  the  case  of  odd  excitation  (Figure  lib),  regions  B* 
and  C  are  excited  in  the  TEM  mode  so  that  the  incident  mode  in  region  B*  is 
out  of  phase  with  the  incident  mode  in  region  C.  Each  of  these  problems  can 
be  solved  individually  for  the  reflection  coefficient  of  the  TEM  mode  in 
region  C.  Let  Re  be  the  voltage  reflection  coefficient  of  the  TEM  mode  in 
region  C  for  the  case  of  even  excitation.  Similarly,  let  Rq  be  the  voltage 
reflection  coefficient  for  the  TEM  mode  in  region  C  for  the  case  of  odd  exci¬ 
tation.  Then  by  super-position, 

"c  ■  \  <Re  r  V  <129> 

R  and  R  can  be  found  quite  simply.  Consider  first  the  case  of  even 
e  o 

excitation.  Refer  to  Figure  12a.  Only  symmetric  transverse  magnetic  modes 
will  be  excited  in  region  A,  i.e.,  the  TM  mode;*  where  n  0  2,  4,  .  „ . 

Because  of  the  symmetry  involved,  it  is  possible  to  place  an  electric  wall  at 
x  -  a/2  and  solve  the  boundary  value  problem  associated  with  the  reduced 
geometry.  This  problem, of  course,  was  solved  in  Section  3.  As  a  matter  of 
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INCIDENT 


^  E  INCIDENT 


REGION  B' 


REGION  B” 


REGION  C 
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x  =  o-c 


x'C 


x  2  0 


Figure  11a.  Even  mode  of  excitation. 
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Figure  lib.  Odd  mode  of  excitation. 
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(a) 

Figure  12a.  Problem  associated  with  even  excitation. 
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XrQ 


(b) 

Figure  12b.  Problem  associated  with  odd  excitation. 
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Then#  after  making  the  appropriate  substitutions  in  the  known  expressions  for 
CC 

S  one  can  write 
oo 


R 

e 


(130) 


Consider  next  the  odd  mode  of  excitation  Pc-fer  to  Figure  12b  Only 
asymmetric  transverse  magnetic  modes  will  be  excited  in  region  Ay  l  e  t  the 
TM  modes  where  n  -  1.  3.  5;  .  - •  A  simpler  but  equivalent  boundary  value 

problem  is  obtained  by  placing  a  magnetic  wall  at  x  -  a/2  The  solution  of 
this  problem  is  quite  straight  forward  In  the  manner  illustrated  in 
Section  3,  a  system  of  infinite  order  linear  algebraic  equations  are  derived 
which  can  be  solved  by  means  of  the  function-theoretic  technique  The  solu¬ 
tion  of  this  problem  yields 


R 

o 


( 132) 
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As  an  exuple  of  the  application  of  Equation  (129),  consider  the  problem 

of  two  equally  spaced  semi-infinite  plates  situation  in  a  waveguide.  Let 

o 

a/X  =  .3.  Using  Equation  (129),  one  computes  Rc  *  .646  e  .  For  ^.rposes 

of  comparison.  Eg  was  also  computed  from  Equation  (127).  With  N  -  4, 

-J49 .0° 

Rc  =  .643e  was  calculated,  which  is  in  very  close  agreement . 

It  should  be  stressed  again  that  Equation  (129)  is  valid  only  for  the 
special  case  of  sysuetrically  spaced  plates.  No  simple  expression  can  be 
found  for  the  more  general  problem  of  arbitrarily  spaced  plates. 
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7  CONCLUSIONS  AND  SUGGESTIONS  FOF  Fl'T'  PF  WORK 
The  generalized  scattering  matrix  technique  has  beer  introduced  and 
applied  to  three  waveguide  discontinuity  problems  They  are  the  E-plane 
metallic  step  discont irui t y,  the  inhomogeneous  E-plane  bifurcation,  and  the 
trifurcated  waveguide..  The  solutions  obtained  arc  formally  exact,  though  in 
series  form. 

i 

V 

To  the  best  of  the  author's  knowledge,  the  general  trifurcation  problem 
has  not  been  discussed  elsewhere.  The  solutions  to  the  other  two  problems, 
viz.,  the  metallic  step  discontinuity  and  inhomogeneous  bifurcated  guide  have 
been  derived,  at  least  approximately,  by  other  authors  using  conventional 
methods  of  analysis. 

Numerical  calculations  have  been  carried  out  for  each  of  the  three  pro¬ 
blems  described  above.  The  results  of  the  numerical  calculations  help  to 

»  ;  '■  , :  , 

demonstrate  the  utility  and  potential  accuracy  of  the  generalized  scattering 

.  .. .  i- 

matrix  technique  applied  to  waveguide  discor tinui ty  problems,  It  is  shown 
that  the  series  expression  for  the  reflection  coefficient;  of  the  dominant 

,  ■  V:,  ■ 

mode  in  the  frequency  range  0  <  a/\  <1,0  is  rapidly  convergent  thiereiby 

making  calculations  convenient..  Moreover,  the  numerical  results.  Where  : 

'•*  i|  ' 

possible,  are  compared  with  published  results  of  other  authors  and  they  are  ; 

••  ■  i| 

found  to  be  in  close  agreement  with  one  another 

The  successful  application  of  the  generalized  scattering  matrix  technique 
to  waveguide  discontinuity  problems  is  a  first  step  in  showing  its  applicabil¬ 
ity  to  a  broader  range  of  problems.  A  list  of  suggested  problems  for  future 
study  together  with  the  corresponding  auxiliary  problems  is  given  below- 

1.  Dielectric  step  discontinuity  in  a  waveguide 

2.  The  diffraction  of  a  plane  wave  by  a  dielect  mc  gi  it  ing  , 


The  suggested  auxiliary  problem  for  1  and  2  is  a  semi -l nf ini te 
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impedance  wall  bifurcating  a  waveguide.  The  impedance  wall  is  char* 
acterized  by  an  impedance  matrix  Z. 

3.  The  diffraction  of  a  plane  wave  by  a  thick  half-plane. 

4.  The  diffraction  of  a  plane  wave  by  a  solid,  circularly  shaped,  metallic 

bar. 

The  auxiliary  problem  is  a  semi-infinite,  tubular  waveguide  in  free  space. 

5.  The  asymmetrical  inductive  and  capacitive  diaphragms  in  a  waveguide  as 

well  as  the  corresponding  strip  grating  problems. 

The  suggested  auxiliary  problem  is  the  bifurcated  waveguide  already  dis¬ 
cussed  in  this  thesis.  Refer  to  Figure  4. 

6.  The  study  of  the  electromagnetic  properties  of  certain  types  of  grating 
structures.  These  structures  have  applications  as  surface  or  leaky  wave 
antennas . 

The  a.tutiliary  problems  for  this  kind  of  problems  is  discussed  in  detail  by 
Ulttra  and  £acewi9 
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